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1.0  Abstract 

Tinder  this  grant  a  study  was  carried  out  of  the  electronic  and 
optical  properties  and  modeling  of  Intercalated  graphite.  Graphite  Inter¬ 
calation  compounds  are  synthetic  metals  with  a  density  much  less  than  that 
of  copper.  By  the  Introduction  of  molecular  Intercalants,  molecular 
conductors  have  been  prepared  which  In  some  cases  have  room  temperature 
electrical  conductivities  comparable  to  that  of  copper.  In  addition  these 
materials  can  have  an  In-plane  to  c-axls  anisotropy  ratio  In  the  electrical 
conductivity  In  excess  of  a  million.  Graphite  intercalation  compounds  are 
therefore  potentially  useful  for  various  applications  in  airborne  systems. 

The  research  program  carried  out  under  this  grant  was  directed  toward  gaining 
an  increased  understanding  of  graphite  intercalation  compounds  In  order  to  be 
able  to  produce  better  materials  with  more  controlled  physical  properties. 

A  major  contribution  of  our  work  has  been  the  formulation  of  a  model  for  the 
electronic  structure,  applicable  to  compounds  of  arbitrary  stage  and  intercalant. 
The  model  has  been  applied  to  the  Interpretation  of  our  Shubnikov-de  Haas 
measurements  which  probe  the  effect  of  Intercalation  on  the  Fermi  surfaces 
of  these  materials.  The  model  has  also  been  extended  to  calculate  for  the 
first  time  the  magnetic  energy  level  structure  of  Intercalated  graphite,  which 
Is  required  to  Interpret  our  magnetoreflection  spectra.  These  spectra  provide 
detailed  and  quantitative  Information  on  the  effect  of  Intercalation  on  the 
electronic  structure  of  intercalated  graphite.  Both  the  magnetoreflection  and 
the  Shubnlkov-de  Haas  experiments  require  high  magnetic  fields  and  utilize  the 
unique  capabilities  of  the  facilities  at  the  Francis  Bitter  National  Magnet 
Laboratory.  Results  on  Infrared  spectroscopy  of  intercalated  graphite  provide 


electronic  structure.  The  use  of  Raman  spectroscopy  to  study  adsorbed 
species  on  graphite  surfaces  Is  demonstrated. 


2.0  Personnel  Involved  with  Research  Program 


M.  S.  Dresselhaus,  Principal  Investigator 
6.  Dresselhaus,  Co-Principal  Investigator 

B. L.  Heflinger,  Graduate  Student  (Surface  Adsorption,  left  September  1979) 
S.Y.  Leung,  Graduate  Student  (Infrared  measurements,  modeling  of 

electronic  dispersion  relations) 

S. A.  Safran,  Post-doctoral  Student  (theory;  left  September  1978) 

C.  Underhill,  Research  Engineer  (experimental  support,  funded  by  MIT 

MauZe  Chair,  left  June  1980) . 

T. C.  Chleu,  Graduate  Student  (Magnetoreflection  experiment,  Landau  level 

calculation) 

M.  Shayegan,  Graduate  Student  (Shubnlkov-de  Haas  experiment  and 
Interpretat ion) 

P,  Haklml  ,  Graduate  Student  (preliminary  Shubnlkov-de  Haas  experiment) 


2.1  Coupling  Activities 

The  MIT  group  is  strongly  coupled  to  international  activities  on  graphite 
intercalation  compounds.  We  are  frequently  invited  to  present  seminars  at 
universities  and  in  industry,  to  give  invited  papers  at  conferences,  and  to 
write  review  articles.  We  collaborate  on  studies  of  the  infrared  electronic 
properties  of  intercalated  graphite  with  Professor  P,C.  Eklund  and  his 
group  at  the  University  of  Kentucky;  they  have  a  strong  experimental  program 
and  we  work  with  them  on  modeling  analysis  of  their  results,  We  have  a 
collaborative  program  with  Professor  J.P.  Issi  at  the  University  of  Louvain- 
la  Hueve  on  thermal  conductivity  and  thermopower  of  intercalated  graphite. 
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with  the  experimental  work  done  by  Issi  and  his  group  and  we  contribute 
to  the  materials  preparation,  modelling  and  interpretation  of  results. 

A  collaborative  program  with  Professor  M.  Corson  and  his  group  at  Boston 
University  has  been  on-groing  on  Mossbauer  studies  of  intercalated  graphite- 
FeCl^  compounds.  We  prepare  their  samples,  they  make  the  Mossbauer  measurements, 
and  we  assist  them  with  the  interpretation  of  their  results. 

2.2  New  Discoveries,  Patents  or  Inventions 
None 


i 
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3.0  Status  of  Research  Effort 

Our  research  program  on  the  electronic  structure  of  graphite  intercalation 
compounds  has  three  major  components.  With  the  magnetoref lection  technique 
we  obtain  experimental  information  highly  relevant  to  the  form  of  the  elec¬ 
tronic  dispersion  relations  and  quantitative  values  for  the  band  parameters. 
With  the  Shubnlkov-de  Haas  techniques  we  obtain  experimental  Imformatlon 
on  the  Fermi  surface.  Information  from  both  of  these  experiments  Is  utilized 
In  our  modelling  calculations  of  the  electronic  structure,  which  Is  valid  for 
any  stage  and  Intercalant.  The  experimental  data  obtained  by  the  Shubnlkov- 
de  Haas  and  magnetoreflection  measurements  are  used  to  determine  the  parameters 
of  the  model  calculations  Including  the  position  of  the  Fermi  level.  Using 
our  model  for  the  dispersion  relations,  the  magnetic  level  structure  (Landau 
levels)  has  been  calculated.  These  Landau  levels  are  now  being  applied 
to  a  self-consistent  Interpretation  of  the  magnetoreflection  and  Shubnlkov- 
de  Haas  experiments.  We  present  below  a  brief  summary  of  our  accomplishments 
on  each  of  these  topics  during  the  grant  period. 

3.1  Studies  of  the  Shubnlkov-de  Haas  Effect  In  Intercalated  Graphite 

Study  of  the  Shubnlkov-de  Haas  (SdH)  effect  In  graphite  Intercalation 
compounds  provides  Information  on  the  Fermi  surfaces  associated  with  each 
of  the  partially  occupied  bands.  During  the  present  grant  period,  a  system 
was  constructed  to  make  Shubnlkov-de  Haas  measurements  on  graphite  Inter¬ 
calation  compounds.  Detailed  measurements  have  been  carried  out  on  a  series 
of  alkali  metal  compounds  with  the  Intercalants  potassium  and  rubidium  and 
these  results  have  been  Interpreted  In  terms  of  a  model  calculation  for  the 
electronic  structure.  This  model  calculation  is  described  in  Section  3.4. 


•MM 
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The  experimental  system  developed  for  carrying  out  these  studies 
is  based  on  the  four-point  method  which  was  set  up  in  a  geometry  which  measures 
the  ii^-plane  transverse  magnetoresistance.  The  system  was  set  up  to  operate 
in  the  temperature  range  1.4<T<4.2 K  and  in  magnetic  fields  up  to  IS  Tesla. 

The  system  was  set  up  so  that  the  angular  dependence  of  the  SdH  oscillations  could 
be  measured  by  rotating  the  sample  around  the  direction  of  the  current  ?  > 
with  I  ^  for  all  rotation  angles,  and  hence  transverse  magnetoresistance  was 
always  measured.  Data  acquisition  was  by  computer  using  a  MINC  system  and  the 
data  were  manipulated  to  obtain  a  Fourier  power  spectrum  of  resistance  vs. 

1/H,  thereby  yielding  directly  the  frequencies  of  SdH  oscillations.  These 
frequencies  are  related  to  the  extremal  cross  sections  of  the  Fermi  surface, 
and  represent  the  final  experimental  output  of  the  measurements. 

During  the  course  of  these  studies,  a  great  deal  of  care  was  given  to 
ensure  the  reproducibility  of  the  data.  For  this  reason,  the  experiment 
was  performed  on  four  potassium  stage  5  samples  and  on  one  of  these,  the 
measurements  were  repeated  several  weeks  later  after  attaching  new  leads. 

The  samples  were  prepared  using  the  two-zone  intercalation  method  with 
HOPG  as  the  host  material.  Because  of  the  instability  of  alkali-metal 
samples  in  the  presence  of  air  and  moisture,  a  great  deal  of  attention 
was  given  to  sample  handling  and  sample  encapsulation.  The  stage  of  the 
samples  was  determined  using  (OOf)  x-ray  diffraction  profiles  both  before 
and  after  the  SdH  experiments,  confirming  that  the  samples  were  single- 
staged  and  Chat  no  desorption  had  occurred  during  the  measurements. 

Our  results  show  that  measurements  made  on  potassium  samples  with 
different  stages  yield  distinctly  different  SdH  frequencies.  Different 
SdH  frequencies  were  also  obtained  with  another  Intercalant,  in  this 


case  rubidium.  Hence  our  conclusion  is  that  there  is  a  unique  set  of 
SdH  frequencies  associated  with  each  stage  and  for  each  donor  Intercalant 
(R,Rb).  Suematsu  et  al.  have  also  shown  stage-dependent  de  Haas-van  Alphen 
frequencies  in  stage  3  and  4  potassium  compounds,  and  these  authors  have 
also  reported  stage  and  Intercalant-dependent  SdH  frequencies.  In  contrast, 
several  authors  have  reported  stage-independent  de  Haas-van  Alphen  type 
frequencies  in  acceptor  compounds,  which  is  inconsistent  with  our  inter¬ 
pretation  of  the  Shubnlkov-de  Haas  results  for  the  alkali  metal  compounds. 

Detailed  results  have  been  obtained  for  graphite- potassium  with  stages  4, 
5  and  8.  These  spectra  have  been  Interpreted  using  our  model  for  the 
electronic  dispersion  relations.  Particular  emphasis  has  been  given  to  the 
interpretation  of  the  results  for  stage  5  graphite-potassium,  where  we  have 
the  most  complete  experimental  results.  Our  model  for  the  electronic 
dispersion  relations  gives  five  valence  and  five  conduction  bands  for  a 
stage  5  compound.  These  bands  are  associated  with  the  graphite  iT-bands, 
and  arise  from  the  k^-axis  zone  folding  procedure.  The  empirical  position 
of  the  Fermi  level  is  determined  to  fit  the  SdH  frequencies,  yielding 
four  partially  occupied  conduction  bands  and  five  completely  filled  valence 
bands.  Various  Fermi  surface  parameters  have  been  calculated  from  this 
model,  including  the  SdH  frequencies,  the  cyclotron  mass  at  E-,  the  electron 

r 

density  for  each  carrier  pocket  and  the  anisotropy  of  the  cross-section 
perpendicular  to  k  .  The  generally  good  agreement  of  the  observed  SdH 

S 

frequencies  with  the  frequencies  calculated  on  the  basis  of  an  empty 
Intercalate  layer  model  suggests  that  the  effect  of  the  intercalant- 
graphite  bounding  layer  interactions  can  be  treated  as  a  perturbation 
and  evaluated  by  fitting  the  model  quantitatively  to  the  observed  SdH 
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frequencles.  Although  the  Shubnikov-de  Haas  frequencies  show  a  significant 
stage  dependence  and  differ  greatly  from  the  SdH  frequencies  for  graphite, 
we  have  shown  that  the  stage  dependence  of  these  frequencies  are  compatible 
with  electronic  dispersion  relations  that  are  basically  graphitic.  The 
application  of  the  band  model  to  explain  the  observed  SdH  frequencies 
has  been  submitted  for  publication  and  a  second  paper  emphasizing  the 
interrelation  between  the  magnetoreflection  and  SdH  results  has  also  been 
submitted  for  publication. 

3.2  Magnetoreflection  Studies  of  Intercalated  Graphite 

We  have  studied  the  electronic  band  structure  of  graphite  intercalation 
compounds  with  the  magnetoreflection  technique  with  particular  reference 
to  spectra  associated  with  interband  Landau  level  transitions.  The 
magnetoreflection  technique  provides  detailed  information  on  the  form  of  the 
electronic  dispersion  relations  and  on  the  values  of  the  band  parameters 
describing  these  dispersion  relations. 

Magnetoreflection  measurements  have  now  been  carried  out  on  well-staged 
acceptor  compounds  (with  the  Intercalants  FeCl^,  AlCl^  and  Br2)  and  on  donor 
compounds  (with  Rb  and  K)  in  a  liquid  helium  cold  finger  dewar  and  in 
magnetic  fields  up  to  15  tesla  at  the  National  Magnet  Laboratory.  Careful 
preparation  of  well-characterized  samples  makes  it  possible  to  satisfy 
the  condition  »  1  and  to  observe  clear  magnetcreflection  oscillations 
for  almost  all  stages  (n  >  3).  The  samples  were  characterized  for  stage 
index  by  (OOL)  x-ray  dlffractograms  both  before  and  after  the  magnetore¬ 
flection  runs.  Typical  spectra  consist  of  a  dominant  structure  that  is 
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identifled  with  Landau  level  interband  transitions  near  the  K-point 
of  the  Brillouin  zone  of  pure  graphite.  A  quantitative  analysis  has 
been  made  of  the  observed  magnetoref lection  spectra,  showing  that  the 
Slonczewski-Weiss-McClure  model  for  pure  graphite  accounts  for  the  dominant 
structure  observed  for  the  graphite  intercalation  compounds,  with  changes 
in  the  band  parameters  that  are  measured  as  a  function  of  intercalant 
and  stage.  Effective  masses  for  the  valence  and  conduction  bands  about 
the  K-point  have  been  determined  and  are  found  to  be  very  small  and  similar 
to  values  obtained  for  pristine  graphite. 

We  found  that  below  a  critical  photon  energy,  which  depends  on  the 
stage  and  the  Intercalant,  no  magnetoref lection  oscillations  are  observed 
experimentally.  This  effect  is  called  the  Fermi  level  cutoff  phenomenon 
and  is  due  to  the  raising  of  the  Fermi  level  in  donor  compounds  and  to  the 
lowering  of  for  acceptor  compounds  as  a  result  of  intercalation.  Thus 
for  the  case  of  donor  compounds,  certain  final  states  will  no  longer  be 
unoccupied  and  for  acceptor  compounds  certain  initial  states  will  no 
longer  be  occupied.  The  Pauli  principle  requires  that  transitions  must 
occur  from  an  occupied  initial  state  to  an  unoccupied  final  state.  Thus 
the  observation  of  the  cut-off  phenomenon  provides  Information  on  the 
position  of  the  Fermi  level  relative  to  the  K-polnt  band  edge  for  valence 
and  conduction  bands  Involved  .with  the  interband  transition.  Thus  the 
cut-off  phenomenon  yields  the  dependence  of  the  Fermi  level  E^  relative 

O' 

to  the  R-polnt  band  edgeg^,^  as  a  function  of  stage  and  Intercalant. 

We  have  obtained  values  for  (Ep-E^,  ^  for  a  number  of  Intercalants  and 
stages. 


The  nagnetoref lection  results  have  been  significant  in  demonstrating 
that  the  electronic  bands  near  the  Fermi  level  remain  highly  graphitic  upon 
Intercalation.  This  conclusion  is  basic  to  our  formulation  of  the  electronic 
dispersion  relations  for  graphite  Intercalation  compounds.  This  work  is 
also  significant  insofar  as  it  yields  quantitative  values  for  the  change 
in  band  parameters  and  Fermi  level  with  stage  and  Intercalant.  These  data 
are  of  Importance  in  obtaining  band  parameters  for  our  electronic  dispersion 
relations.  Our  present  efforts  focus  on  the  calculation  of  the  magnetic  energy 
level  structure  (Landau  levels)  using  our  model  for  the  electronic 
dispersion  relations.  These  Landau  levels  will  then  be  supplied  to  yield 
a  consistent  Interpretation  of  the  observed  magnetoreflection  spectra,  since 
the  present  analysis  of  the  magnetoreflection  results  has  been  on  the 
basis  of  a  graphite  Tr-band  model. 

Since  the  magnetoreflection  results  for  intercalated  graphite  are 
very  similar  to  those  for  pristine  graphite  while  the  Fermi  surface 
measurements  show  that  major  changes  occur  in • the  Fermi  surface  upon  inter¬ 
calation,  it  was  believed  by  some  authors  that  the  magnetoreflection  and 
Shubnlkov-de  Haas  results  are  inconsistent.  Our  work  however  shows  that  the 
results  from  both  experiments  can  be  explained  on  the  basis  of  energy 
bands  that  are  highly  graphitic,  but  that  energy  level  shifts  and  Fermi  level 
shifts  occur  as  a  result  of  Intercalation.  A  paper  showing  the  consistency 
of  the  magnetoref lection  and  Shubnlkov— de  Haas  results  has  been  prepared  and 
will  appear  in  the  Proceedings  of  the  15th  International  Conference  on 
the  Physics  of  Semiconductors.  An  invited  talk  summarizing  the  magnetore¬ 
flection,  Shubnlkov-de  Haas  results  and  their  interpretation  based  on  our 
band  model  has  been  presented  at  Conference  on  the  Application  of  High 
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Magnetic  Fields  Co  Che  Physics  of  Semiconduccors  ac  Hakone,  Sepcember 
1980,  and  will  appear  in  Che  Conference  Proceedings,  published  by  Springer. 
We  are  che  only  group  Incernaclonally  chac  has  presenced  resulcs  on  the 
raagnetorefleccion  spectra  of  graphite  intercalation  compounds. 

3.3  Calculations  of  the  Electronic  Energy  Band  Structure  of  Graphite 

Intercalation  Compounds 

We  have  calculated  che  electronic  energy  bands  for  graphite  inter¬ 
calation  compounds  using  a  model  based  in  che  well-established  electronic 
structure  of  pristine  graphite.  In  our  work,  we  used  che  Hamiltonian 
developed  by  Slonczewski-Weiss-McClure  (SWMcC)  for  che  graphite  ir-bands, 
and  have  folded  these  bands  along  the  k^  axis  to  account  for  the  super- 
lattice  periodicity  introduced  by  staging,  The  formulation  can  also  treat 
an  in-plane  superlatcice  periodicity  by  In-plane  zone  folding,  in  cases 
where  such  periodicity  is  found  experimentally.  For  convenience  in  intro¬ 
ducing  che  perturbation  due  to  the  intercalant,  che  folded  k-space  Hamil¬ 
tonian  is  transformed  into  a  layer  representation.  To  obtain  the  proper 
staging  periodicity  we  replace  one  layer  out  of  every  (n+l)  graphite 
layers  by  an  Intercalant  layer.  With  the  inclusion  of  an  interaction 
between  the  intercalate  layer  and  the  graphite  bounding  layer,  electronic 
structures  appropriate  to  both  donor  and  acceptor  compounds  can  be  obtained. 
The  model  is  applicable  to  any  intercalant  and  any  stage,  and  is  therefore 
ideally  suited  to  che  interpretation  of  many  experiments.  In  our  research 
program  we  have  explicitly  applied  the  model  to  the  interpretation  of 
che  magnetoref lection  experiment,  which  is  sensitive  to  the  magnetic 
energy  levels  and  che  electronic  dispersion  relations,  and  to  che  Shubnikov- 
de  Haas  experiment  which  is  sensitive  to  che  Fermi  surface.  Reflectivity 
experiments  are  planned  for  future  work  to  yield  structure  associated  with 
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electronlc  transitions  that  can  be  used  to  test  the  model  developed  for  the 
electronic  structure  on  the  basis  of  magnetoref lection  and  Fermi  surface 
measurements . 

One  paper  on  the  general  formulation  of  this  model  is  now  in  press. 

A  second  paper  which  applies  this  model  to  the  interpretation  of  Shubnikov- 
de  Haas  results  was  presented  at  the  Provincetown  Conference  and  will  be 
published  in  Synthetic  Metals.  A  third  paper  which  contains  details  on  the 
treatment  of  the  In-plane  superlattice  symmetry  was  presented  at  the 
Yamada  Conference  on  Layered  Materials  and  will  be  published  in  Physica. 

3.4  Calculation  of  the  Magnetic  Energy  Level  Structure  for  Graphite 
Intercalation  Compounds 

Using  the  basic  Hamiltonian  for  the  electronic  dispersion  relations 

described  in  Section  3.3,  we  have  for  the  first  time  calculated  the 

magnetic  energy  level  structure  for  a  graphite  intercalation  compound. 

In  this  calculation  the  Luttinger  transcription  for  the  wave  vector  operator 

k  was  used,  and  only  linear  terms  in  k  were  retained  in  the  off-diagonal 

terns.  Trigonal  warping  was  also  neglected  to  slraplifly  the  calculation. 

With  these  approximations  the  Landau  levels  can  be  obtained  in  closed  form. 

Explicit  results  have  so  far  been  obtained  for  stage  1,3  and  5  compounds. 

Zone  folding  along  k  causes  a  near  degeneracy  of  a  number  of  bands  in  the 

z 

vicinity  of  the  HK  axis,  and  this  gives  rise  to  crossings  of  the  magnetic 
energy  levels.  Results  are  now  being  obtained  for  higher  stage  compounds 
for  which  we  have  detailed  magnetoreflection  results  so  that  the  calcu¬ 
lations  can  be  applied  directly  to  the  interpretation  of  the  experimental 
data.  Perturbations  due  to  higher  terms  in  k,  trigonal  warping  and  inter¬ 
action  between  the  intercalant  and  the  graphite  bounding  layer  will  be 
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included.  This  will  enable  us  to  obtain  Information  on  the  magnitude  of 
these  terms,  which  will  be  utilized  In  the  model  for  the  electronic 
structure. 

3.5  Preparation  of  Review  Article 

During  the  period  of  this  grant,  a  monograph  was  prepared  on  the  general 
field  of  graphite  intercalation  compounds.  The  largest  chapter  of  this  monograph 
deals  with  the  electronic  properties,  and  was  supported  by  this  grant.  This 
monograph.  Invited  by  Advances  in  Physics,  has  been  completed  and  Is  now  In  the 
review  process. 

3.6  Infrared  Studies  of  Graphite  Intercalation  Compounds 

During  the  period  July  15,  1978  -  July  14,  1979,  this  grant  supported  in¬ 
frared  studies  on  intercalated  graphite.  During  this  period  we  studied  the 
infrared  spectra  of  graphite  Intercalated  with  the  acceptors  FeCl^  and  AlCl^ 
and  with  the  alkali  metal  donors  K,  Rb  and  Cs.  The  spectra  were  taken  on 
a  Fourier  Transform  Spectrometer,  provided  by  a  Central  Facility  of  the  Center 
for  Materials  Science  and  Engineering. 

Our  Infrared  studies  of  the  lattice  modes  of  graphite  intercalation  com¬ 
pounds  have  shown  that  the  odd  parity  of  the  graphite  modes  is  preserved  upon 
intercalation  and  that  the  oscillator  strength  provides  significant  information 
on  the  dynamic  effective  charge  and  the  magnitude  of  the  dipole  moments 
associated  with  the  lattice  mode.  From  analysis  of  the  dependence  of  the 
IR-active  mode  frequencies  on  reciprocal  stage,  we  have  been  able  to  identify 
modes  that  are  associated  entirely  with  graphite  interior  layers  and  modes 
which  involve  graphite  bounding  layers.  For  the  donor  compounds  the 
infrared-active  mode  frequencies  characteristically  decrease  as  a  function 
of  reciprocal  stage  (1/n)  while  for  acceptor  compounds,  the  mode  frequencies 
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increase  as  a  function  of  (1/n).  A  similar  behavior  is  found  for  the  Raman 
active  modes,  consistent  with  the  observed  expansion  of  the  In-plane  lattice 
constant  for  donor  compounds  as  a  function  of  (1/n) ,  and  suggesting  that  the 
in-plane  lattice  constant  decreases  as  a  function  of  (1/n)  for  acceptors.  This 
prediction  has  now  been  confirmed  experimentally  in  work  on  graphite-NiCI^ 
compounds  by  Flandrols  and  coworkers  at  Bordeaux  and  also  in  work  in  our  own 
group  by  Tania  Krapchev  on  graphlte-FeCl^. 

Our  infrared  studies  have  focussed  on  the  similarities  and  differences 
in  behavior  of  the  infrared  lattice  mode  spectra  for  donors  and  acceptors. 

We  have  found  that  infrared-active  modes  are  not  observed  for  either  donor 
or  acceptor  compounds,  as  expected  since  two  layers  are  required  to  produce 
an  Infrared-active  mode.  For  stage  2  compounds,  the  acceptors  are  infrared- 
active,  but  the  donors  are  not.  The  stage-dependent  frequency  shift  as  a 
function  of  reciprocal  stage  is  positive  for  the  acceptors  and  negative  for  the 
donors.  The  magnitude  of  the  frequency  difference  between  modes  involving 
only  interior  layers  and  modes  Involving  bounding  layers  is  much  larger  for 
donor  compounds  than  for  acceptors.  The  relative  oscillator  strengths  of 
interior  and  bounding  layer  modes  indicate  that  more  charge  is  transferred 
to  the  graphite  interior  layers  for  the  donors  than  for  the  acceptors, 
consistent  with  electrical  conductivity  and  magnetoref lection  results. 


3.7  Study  of  adsorbed  Molecules  on  Graphite  Surfaces  Using  Raman  Spectroscopy 

This  study  was  initiated  to  study  the  connection  between  surface  adsorption 
and  intercalation.  Bromine  was  chosen  for  the  adsorbing  species  because  Br^ 
is  the  only  molecular  intercalant  for  which  Raman-active  in-plane  modes  have 
been  observed  and  because  the  spectra  of  the  free  molecule  have  been  studied 
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extensively. 

The  Raman  spectra  at  room  temperature  consist  of  numerous  sharp  lines 
in  the  wavenumber  region  290  <u)< 340cm  Extremely  strong  resonant  Raman 
effects  were  observed  when  the  frequency  of  the  incident  light  was  varied  by 
means  of  a  tunable  dye  laser .  A  given  Raman  spectrum  changes  radically  by 
varying  the  laser  energy  by  as  little  as  1  cm~^.  These  complex  spectra  can 
be  understood  by  a  resonant  Raman  process  involving  vibrational-rotational 
levels  of  a  free  Br2  molecule  in  the  ground  or  first  excited  electronic 
state. 

The  Raman  spectra  observed  on  a  low  temperature  99<T<230K  graphite 
surface  were  distinctly  different  from  the  resonantly-enhanced  free  molecular 
spectrum.  The  low  temperature  spectra  showed  b>ro  or  more  very  broad  lines 
distinctly  different  from  either  those  for  the  free  molecule  or  for  the 
Intercalated  Br^  molecule.  The  low  temperature  Raman  spectra  were  thus 
attributed  to  adsorbed  bromine  species.  The  Ph.  D.  thesis  of  Bruce 
Hef linger  was  based  on  this  work. 
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Raman  Study  of  Adsorbed  Br?  on  fraphlte.*  B . L . 
MEFLINGER,  M.S.  DRESSELHAUS  and  u.  DRESSELRXUSt,  MIT— 
Bromine  intercalates  into  graphite  abov«»  a  threshold 
vapor  pressure  of  18  Torr  at  rscra  temp.^  At  lover 
pressure  bromine  uptake  is  by  surface  adsorption.  We 
report  Raman  spectra  from  high  surface  area  forms  of 
graphite  (Crafoil  and  exfoliated  graphite)  at  vapor 
pressures  ranging  between  10~3  Torr  -«nd  20  Torr,  where 
the  transition  from  surface  adsorption  co  balk  inter¬ 
calation  occurs.  We  find  the  Raman  spectra  co  be 
sensitive  to  both  surface  and  bulk'lattice  modes.  The 
surface-specific  modes  are  shifted  60  cm“i  higlier  in 
frequency  than  the  bulk  graphite  modes  (1580  cm“l  and 
1600  cm“l)  found  in  the  intercalation  compounds.  We 
.also  report  the  frequency  dependence  of  the  bromine 
stretch  mode  as  a  function  of  surface  coverage. 

*  Sponsored  under  AFOSR  Contract  f 77-3391 
t  Natl.  Magnet  Lab,  supported  by  NSF 
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Prefer  Regular  Session  Submitted  by 

•  « 

?yi ,  ^  /Pi:  3  o«-  t.A  u 

M.S.  Drcssclhnus 
Room  13-2090,  M.I.T. 

I  Cambridge,  M\  02139 

Note:  To  be  presented  after  p.npcr  by  P.C,  Eklund,  N.  Kambe,  G.  Dresselhnus 
and  M.S.  Drcssclhaiis  entitled  "Resonant  Raman  Scattering  From  nr2 
in  Graphite  Intcrc.il.'ttlon  Compounds". 


Abstract  Siibir.lLCed 


/ 


for  the  Chicago  Kucclng  of  the 


'Physical  Reviev; 
Classification  Scheme 
rHu^ber  78.20 


American  Physical  Society 


Bulletin  Subject  Head  ini; 
Graphite  Tntcrc.ilnti on 
Compounds 


Fermi  Surf.ace  Model  for  Graphite  FeClS.t  C. 
DRESSELIL\US*and  S.Y.  LEUNG,  M.I.T.--The  angular  dependence 
of  the  various  Shuhnikov-de  Haas  (SdH)  periods  on  staRC 
Is  interpreted  in  terms  of  fluted  graphitic  ellipsoidal 
Fermi  surfaces  which  have  extremal  cross-sectional  areas 
.resulting  from  the  staging  of  the  intercalation  compounds. 
The  structure  factor  calculations  indicate,  in  agreement 
with  the  interpretation  of  0-  29  x-ray  spectra,  that 
although  the  folded  Erillouin  zone  for  the  dilute  compouds 
is  very  small  (indicating  many  extremal  areas)  only  the 
Bragg  reflections  near  those  for  pristine  graphite  have 
structure  factor  sufficiently  large  to  yield  observable 
SdH  periods.  This  model  accounts  for  (1)  the  observed 
•  angular  dependence  and  cut-off  phenomena  for  the  slow  a>id 
fast  periods,  (2)  the  weak  stage  dependence  of  the  fast 
periods  and  (3)  the  larger  stage  dependence  of  the  slow 
periods.  The  model  also  explains  why  the  anisotropy  of 
the  Fermi  surface  is  independent  of  stage  and  of  the 
volume  enclosed  by  the  constant  energy  surfaces  and  is 
similar  to  that  for  pristine  graphite. 


t  Supported  by  AFOSR  Grant  f77-3391. 


*  Kat'l  Magnet  Lab,  supported  by  NSF. 
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Infrared  Spectrosconv  of  Graph ire-Alkali  Metal 
Intercalation  Coffno>nd.V'4:  S.V.  LEUICG,  C.  L’NDEiUnLL,  G. 
DRESSELllAUfiv  M .  I . T . — Tor  the  first  time,  infrared  spec¬ 
tra  of  graphite-alkali  metal  donor  compounds  are  repor¬ 
ted.  Of  significance  is  the  large  ('  30cn~i)  do'.'nshift 
of  the  infrared-active  optic  modes  from  the  r.]^u  mode  in 
graphite  at  1537cia"i,  in  censtrast  with  results  observed 
for  the  acceptor  compounds  v;here  the  stage  dependent 
frequency  shift  is  smaller  in  magnitude  and  opposite  in 
sign.  Results  are  reported  for  the  IR-active  mode 
frequencies,  linewidths  and  oscillator  strengths  for  K, 
Rb,  and  Cs  compounds  as  a  function  of  stage.  Analysis 
of  the  oscillator  strengths  nornits  independent  determi¬ 
nation  of  the  charge  transfer  to  the  interior  and  ' 
bounding  graphite  layer  modes  vs.  stage  for  K,  Rb  and 
Cs  donor  compounds.  Analysis  of  the  stage  dependence  of 
the  mode  frequencies  indicates  an  in-planc  lattice 
expansion  due  to  the  intercalation  process,  in  agreement 
with  direct  x-ray  measurements.* 

•d.E.  Nixon  and  G.S.  Parry,  J.  Phvs.  C.  2,  1732  (1969). 
tSupported  by  AFOSR  Grant  ^-77-3391. 

*Nat'l  Magnet  Lab,  supported  by  NSF. 
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Direct  Determination  of  Charge  Transfer  In 
Graphite  Intercalation  Compounds."  M.S.  DRESSCLHAUS, 
■S.Y.  LEUiNG,  C.  UNDERHILL  and  G.  DRESSELIL'lUS'*’,  MIT  — 
Although  charge  transfer  in  graphite  acceptor  compounds 
has  been  discussed  by  many  authors,  only  recently,  using 
infrared  spectroscopy,  has  the  charge  on  the  graphite 
bounding  and  interior  layers  been  determined  indepen¬ 
dently.  This  determination  is  possible  because  the- IR 
modes  associated  with  the  interior  and  bounding  layers 
occu]^  at  different  frequencies,  and  a  lineshape  analysis 
of  the  IR  spectra  yields  the  dipole  moment  and  hence 
charge,  associated  with  each  mode.  For  example,  the 
results  for  graphite-FeClj  show  that  there  is  a  70% 
charge  transfer  to  the  bounding  layers  (essentially 
independent  of  stage),  and  for  stage  n>.3  an  additional 
30Z  transfer  is  made  to  the  interior  layers  collectively. 
The  small  charge  transfer  to  graphite  interior  layers  is 
corroborated  by  magnetoreflection  results  on  graphite- 
FeCl3  which  probe  primarily  the  graphite  interior  layers. 
We  also  discuss  the  relation  of  this  determination  of 
the  charge  transfer  to  Shubnikov-de  Haas  data,  which 
yield  fluted  Fermi  surface  ellipsoids  associated  with 
both  graphite  interior  and  graphite  bounding  layers. 
tNat'l  Magnet  Lab,  supported  by  NSF. 
i*Supported  by  AFOSR  Grant  #77-3391. 
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I  Infrared  and  Raman  active  nodes  in  graphite 
intercalation  conpounds  can  be  classified  into  j 
fchree  categories:  those  arising  troa  (1)  the  ; 

intercalant;  (2)  the  graphite  bounding  layers  (the 
graphite  layer  that  is  adjacent  to  the  Ictcrcalant) 
and  (3)  the  interior  graphite  layers. Since  the 
^ttlee  nodes  of  the  Intcrcalanc,  graphite  bounding 
had  interior  layers  occur  at  different  frequencies, 
these  techniques  can  investigate  each  of  the 
different  types  of  layers  separately.  The  present 
study  focusses  on  two  Important  Ideas  in  the  under¬ 
standing  of  graphite  Intercalation  conpounds: 

Stijln  field  and  effective  charge.  We  have  ' 

■ensured  the  stage  dependence  of  the  frequency 
Upshift  of  the  graphitic  optic  E2g2  ^lu  nodes 
l(■rialDg  from  both  the  bounding  and  interior 
graphite  layers).  A  sinilar  stage  dependent  fre¬ 
quency  upshift  Is  observed  In  graphite  ferric 
chloride  for  all  the  infrarcd-accivc  and  Raaan- 
■ctlve  lattice  nodes.  This  stage-dependent  upshift 
bf  the  lattice  codes  can  be  interpreted  via  a 
Strain  field  Intercalation  r.cchanisn.  Analysis  of 
the  Infrared  spectra  gives  the  stage  dependence  of 
the  effective  charge  associated  with  the  graphitic 
founding  and  interior  layers.  An  effective  charge 
rstln  of  7/3  was  measured  between  the  bounding  and 
'  Interior  layers  for  graphite-ferric  chloride. 

I  Vslog  highly  oriented  pyrolytic  graphite  (ItOPC) 
as  s  host  naterlal  and  employing  the  conventional 
two-zone  technique^  we  have  been  able  to  grow  a 
series  of  essentially  single-staged  graphite-ferric 
chloride  samples  including  stages:  1,2, 3, A, 5, 6,  and 
11.  The  stage  of  the  corqtound  was  determined  using 
S  conventional  0-29  x-ray  diffractometer  scan 
caploying  MoKq  as  the  incident  radiation.  Adolx- 
tures  of  secondary  stage  material  appear  as 
additional  peaks  In  the  x-ray  spectra  and  these 
additional  peaks  can  easily  be  Identified.  | 

!  I 

I  Koom  temperature  Raman  spectra  were  t.nken 
using  an  argon-ion  laser  operating  at  A880A  in  the 
bacfc-acatterlng  geometry.  With  the  incident  beam 
along  the  c-axls  of  the  crystal,  the  Raman-active  i 
in-plane  E2g2  modes  are  excited.  Intcrcalant 
desorption  associated  with  laser  heating  was  > 

avoided  by  using  a  low  laser  power  level  «  SO  nw). 
Infrared  reflectivity  spectra  were  taken  using  a 
lourler  transform  spectrometer.  Room  temperature 
apeetra  covering  the  energy  range  AOOfuFAOOO  cm~l 
were  taken  with  the  sample  e-axls  parallel  to  the 
Incident  beam.  I 


I  Raaan  and  Infrared  spectra  of  graphltc-ferrlc  - 
dilorlde  are  consistent  with  the  Identification  of 
the  lattice  modes  with  'he  graphite  bounding  and 
Interior  layers. ^  In  the  Raman  spectra  a  tingle 
borcntzlan  line  is  found  for  both  stage  1  and  2 
ceapounds  and  two  Lorcntzlan  lines  for  compounds 
ylth  n>l.  For  the  Infrared  spectra,  fitting  the  data 
!■  tens  of  a  Drude  background  and  Lorcntzlan  oscil¬ 
lators  also  yields  a  alngle  Lorcntzlan  oscillator 


for  stages  1  and  2  and  two  Lorcntzlan  oscillators 
•for  stages  n^l.  The  higher  frequency  component 
In  the  Raman  spectrum  (E2g2)  is  identified  with 
the  bounding  graphitic  layer  because  It  exists  in 
all  stages  (including  1  and  2)  but  not  In  pristine 
graphite,  while  the  lower  frequency  component 
(E2g,)  is  Identified  with  the  interior  graphite 
layer  because  its  frequency  extrapolates  to  that  for 
pristine  graphite  in  the  Unit  of  very  diluC’i 
compounds.  Furthermore,  there  is  no  evidence  of 
the  lower  frequency  component  in  stages  1  and  2. 
Identification  of  the  infrared-active  rodes  with 
the  bounding  and  interior  layers  follo-ws  a  similar 
line  of  reasoning.  The  difference  is  that  for  the 
Infrared  spectra  the  high  frequency  component  (t-.-j) 
is  identified  with  the  interior  layer  and  the  It-.cr 
frequency  component  (Elu)  is  identified  with  the 
hounding  graphite  layer. 

The  dependence  of  the  lattice  node  frequency  u 
on  reciprocal  stage  (1/n)  for  both  the  infrared- 
active  and  Raman-active  modes  is  shewn  in  Fig.  1. 

We  observe  that  all  the  lattice  node  frequencies 
(infrared  and  Raman)  for  n%3  exhibit  the  same 
linear  dependence  on  reciprocal  stage  whether  the 
mode  is  associated  with  graphite  bounding  or 
interior  layer. 

The  frequency  upshift  of  the  modes  shown  in 
Fig.  1  is  Interpreted  as  a  lattice  stiffening 
upon  Intercalation.  As  the  number  of  graphite 
layers  between  successive  intercalate  layers 
Increases,  the  strain  experienced  per  graphite 
layer  Is  decreased' and  thus  a  smaller  shift  in  the 
lattice  mode  frequency  results.  The  slope  of  the 
dependence  of  the  mode  frequency  on  reciprocal 
■Cage  has  the  seme  sign  in  graphite-AiCi3^®,  a 
compound  similar  in  structure  to  graphite-FeCij 
(l.e.  for  both  types  of  acceptor  conpounds  the 
mode  frequencies  increase  as  a  function  of  1/n). 

In  contrast,  for  both  the  infrared-active  modes*! 
and  Raman-Bctlve3  modes  in  the  graphite-alkali 
metal  conpounds,  the  slope  of  a  plot  of  the  mode 
frequency  vs  reciprocal  stage  has  a  sign  opposite 
to  that  for  graphltc-ferrlc  chloride  (l.e.  the 
■ode  frequencies  decrease  as  a  function  of  1/n). 
This  Interpretation  of  the  mode  frequency  shift  Is 
consistent .with  the  observation  that  in  graphite- 
potassium  compounds  there  Is  a  contraction  in  the 
graphite  In-plane  lattice  constant  (ao)  upon 

Intercalation. i 


I  Infrared  spectroscopy  la  useful  for  the  study 
of  lattice  nodes  and  electronic  properties  of 
graphite  Intercalation  compounds  in  the  following 
two  ways:  (1)  infrared  spectroscopy  probes  o. d- 
parlty  nedee  while  Ranan  spectroscopy  probes  cven- 
parlty  nodes  In  these  compounds.  (2)  Detetnination 
of  the  oscillator  strengths  provides  Important 
Information  concerning  electronic  properties  such 
as  the  dipole  moments  and  effective  charges.  »he 
dipole  nonent  associated  with  the  Inf rarcd-aetlvc 


^Ihe  Raman  and  Infrared  studies  were  respectively  supported  by  the  Office  of  Naval  Research  (ONR  .V00014-77- 
*003)  and  the  Air  Force  Office  of  Scientific  Research  (f 77-3391). 

fpepertoent  of  Electrical  Engineering  and  Computer  Science  and  Center  for  Meterlala  Science  and  Engineering, 
irr*««**  litter  National  Magnet  Laboratory,  supported  by  the  NSF.  - — - — -  '  ' 


ao  a  {unction  o(  reciprocal  s'.ar.n  is  shovn  I 
as  en  tnnot  to  fip,.  1  for  bouniiinr  mid  Interior 
(raplilte  layers,  blnce  difCccunt  types  of  graphite 
layers  can  be  probed  Indupcndcntly,  ve  have  | 

datetstned  the  effective  chatgc  associated  ulth  the 
koundlnc  and  interior  layers  separately.  Analysis 
of  the  dipole  eonent  in  terms  of  an  effective 
eharse^  shows  that  for  the  graphice-FeCl^  systea 
the  charge  distribution  is  divided  such  that  702  of 
the  charge  resides  on  the  two  bounding  layers  and 
the  reaalning  302  is  shared  among  the  n-2  interior 
layers.  These  charge  distributions  will  be 
related  to  (1)  Fetal  surface  measurements  on  the 
seas  compounds  using  the  Shubnikov-de  Haas  Efface 
(reported  by  Voollaa  ec  al  at  this  conference) 

(2)  the  Femi  level  shift  in  the  graphite  interior 
layers  as  a  result  of  Intercalation  with  FCCI3  as 
Measured  by  the  nagnctoreflcctlon  technique 
(reported  by  Mendct  ee  al  at  this  conference)  and 

(3)  the  dependence  of  the  electricel  conductivity 

ea  reciprocal  etage.  | 
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.  Pig.  1  The  inverse  stage  dependence^of  the 
node  frequencies  for  Raman  Ejgj  ^<*4  E2g2 
(squares)  and  IR  Exu  <'>hd  £xu  (circles)  modes 
In  graphite-ferric  chloride  compounds.  The  * 
and  superscripts  refer  respectively  to 
graphite  Interior  and  bounding  layers.  The 
solid  lines  represent  a  lease  squares  fit  to 
the  experimental  points  for  n%3  and  Indicate 
S  similar  strain  dependence  for  all  observed 
•odes.  The  Ejgj  And  Ex^  mode  frequencies 
lor  pristine  graphite  arc  given  at  (l/n)>0. 
Ihs  inset  shows  c.-te  (1/n)  dependence  of. the 
dipole  moment  (or  the  Infrared  Eiu  and  Exu 
■odes  noraallzcd  to  chat  for  the  Eiu  mode  In 
pristine  graphite.  Tlic  Raman  and  IR  frequen- 
slss.lor  the  lewesj  stages  q>l  and  n>2 
compounds  labeled  E2gs  *"4  Exu  depart  from 
the  hahavler  obaacved'for  the  niJ  compounda. 
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Recent  experlocntal  results  on  the  electronic 
atructure  of  graphite  intercalation  compounds  such 
as  those  reported  at  this  conference  using  the 
Shubnikov-de  Haas  effect^  and  the  magneto-  . 
reflection  technique  2  show  that  the  electronic 
levels  near  the  Fermi  level  are  mainly  derived 
from  graphite  s-bands.  These  experimental 
programs  are  currently  generating  large  amounts  of 
experimental  data  relevant  to  the  dependence  of 
'the  electronic  structure  on  intercalate  species 
land  stage,  and  there  is  consequently  a  pressing 

an  be  used  to 
.  First 
he  stage  1 

jalkali  metal  compounds  CgLi^  and  CgK^,  and  these 
■calculations  confirm  that  the  significant  energy 
'bands  near  the  Fermi  level  originate  from  the 
'graphite  x-bands.  On  Che  ocher  hand,  an  excensive 
amount  of  work  is  required  to  determine  EOt)  from 
!a  first  principles  calculation  for  each  intercalate 
species  and  stage,  and  the  calculations  become  more 
difficult  with  increasing  stage.  Because  of  the 
success  of  phenomenological  models  based  on 
.aymmetry  (e.g.  the  Slonczewski-Ueiss-McClure 
(S-W-McC)  model)^,b,7  Che  interpretation  of 
■experimental  data  for  graphite,  it  would  seem  that 
’Useful  phenomenological  models  based  on  the 
jSyanietry  of  the  Intercalation  compounds  could  also 
play  an  Important  role  in  the  Interpretation  of 
■experimental  data  for  these  materials. 

I 

!  Vc  present  here  a  phenomenological  model  for 
'the  calculation  of  the  electronic  dispersion 
.relations  for  graphite  intercalation  compounds 
baaed  on  the  symnetry  properties  of  the  graphite 
structure  upon  which  the  superlattlcc  symmetry  of 
;the  Intercalation  compounds  is  superimposed.  The 
model  Is  readily  applicable  to  intercalated 
■graphite  for  any  intercalate  species  and  stage. 

;la  the  limit  of  the  stage  1  alkali  metal  compounds 
CgLl  and  CgK,  the  dispersion  relations  calculated 
■on  the  basis  of  the  phenomenological  model  are  in 
'good  agreement  with  those  reported  using  a  first- 
principles  calculation.  In  the  limit  of  very  high 
Stage,  the  S-W-McC  graphite  model  is  recovered  for 
jk  points  near  the  edges  KKH  of  the  Brlllouin  zone. 
The  basic  features  of  this  band  structure  model  are 
determined  directly  from  the  band  parameters  of 
pristine  graphite  (which  are  known)  and  the 
auperlattlce  periodicity  (which  is  deduced  from  x- 
ray  data).  A  quantitative  fit  of  the  model  to  the 
experimental  data  now  emerging  from  the  Fermi 
,s(irfacc  and  magnetoref lection  spectra  will 
determine  the  band  parameters  representing  Inter- 
jacclon  between  the  graphite  and  intercalate  layers. 

I  The  phenomenological  Mdcl  presented  here  tukes 
mao  of  the  3-dlmcnsional  Fourier  expansion  of  the 
grephlte  x-bands  previously  developed  by  Johnson 
mod  Drcsselhaus^  to  account  for  the  optical 
properties  of  pristine  graphlec.  This  model  Is  an 
I 


!need  for  energy  band  models  which  c 
llnterpret  these  data  quantitatively 
iprlnclples  calculations  exist  for  t 


extension  of  the  S-V-McC  model  which  Is  based  on  a 
Foj^r^cr  expansion  of  the  energy  bands  along  kz  and 
a  k>p  expansion  about  the  HKH  axes.  Thus,  by 
carrying  out  the  Fourier  expansion  in  3-dincnsions , 
a  model  Is  obtained  for  E(k)  throughout  the 
Brlllouin  zone.  This  extension  is  necessary  for 
the  interpretation  of  the  optical  data  which 
involves  contributions  from  a  larger  volume  of  the 
Brlllouin  zone  than  does  Feral  surface  data.  To 
carry  out  the  zone  folding  procedures  outlined 
below,  such  an  extension  is  also  necessary'  for 
application  to  the  intercalation  compounds.  It 
should  be  emphasized  chat  the  use  of  the  3- 
dimenslonal  Fourier  expansion  for  graphite  does  not 
introduce  any  parameters  other  than  those  already 
evaluated  'oy  Johnson  and  Dressclhaus  in  their 
Interpretation  of  the  optical  data  for  graphite.^ 

The  present  phenomenological  model  adds  to  the 
basic  graphite  symmetry  the  superlattice 
periodicity  appropriate  to  the  intercalation 
compounds  through  suitable  zone  folding  of  Che 
basic  graphite  Brlllouin  zone.  The  in-plane 
intercalate  order  gives  rise  to  zone  folding  of 
the  in-plane  wave  vectors  and  in  this  way  is 
sensitive  to  the  intercalate  species.  The  staging 
phenomenon  introduces  a  c-axls  superlatcice  struc¬ 
ture  which  is  incorporated  into  the  calculation  by 
C-axis  zone-folding.  The  zone  folded  representa¬ 
tion  is  then  transformed  into  a  layer  representa¬ 
tion  and  the  interaction  between  graphite  bounding 
layers  and  the  intercalate  layer  Is  treated  as  a 
perturbation.  This  perturbation  is  dependent  on 
band  parameters  that  are  related  to  the  overlap 
between  the  graphite  n-band  orbitals  and 
orbitals  on  the  intercalate  layer.  These  band 
parameters  can  be  evaluated  by  comparison  with 
either  first  principles  calculations  or  experimen¬ 
tal  Fermi  surface,  optical  and  magneto-optical 
data.  It  should  be  noted  Chat  once  these  inter¬ 
action  band  parameters  are  evaluated  for  a  specific 
Intercalant,  the  phenomenological  bands  for  that 
Intcrcalant  are  determined  for  all  stages  that 
have  the  tame  in-plane  ordering. 

The  results  for  E(k)  for  the  phenomenological 
model  for  stage  1  compounds  with  the  CgX  and  CgX 
are  shown  in  Figs.  1  and  2  using  no  adjustable 
parameters.  In  this  form,  the  bands  are  broadly 
applicable  Co  all  compounds  with  the  indicated 
structures.  For  each  figure  the  results  are 
compared  with  the  pertinent  first-principles 
calculations  for  CgLl}  and  CgX^  and  very  good 
agreement,  is  obtained  for  Che  bands  related  lj 
the  graphite  n-bandi.  To  obtain  agreement  for  the 
hybfldlzcd  metal  s  and  p  bands,  it  is  necessary  to 
fit  the  phenomenological  bands  to  the  first 
principles  hybridized  s-bands  to  determine  the 
overlap  band  parameters  between  the  graphite  p, 
w-orbltals  and  the  intercalate  orbitals.  Ono 
attractive  feature  of  the  phenomenological  model  In 
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that  once  the  band  pnraneCcra  arc  dpcerslncd  for 
any  low  sense  conpoimd  (where  the  (’.rnphlte* 
Intercalate  Interact  tons  are  Inportnnt),  the 
dlsper  ston  relations  for  all  stases  can  then  be 
evaluated  without  introduction  of  additional  band 
parancters.  The  phonomcnolostcal  nodcl  thus 
provides  a  useful  Cool  for  Ct>e  calculation  of  the 
electronic  structure  for  high  senp,c  compounds 
where  first  principles  calculations  are  difficult 
to  Bake  and  the  phenomenological  model  Is  expected 
to  be  Bost  convergent. 


Ce  X  Stage  1 

ALMTKHAr 


Fig.  1  Coaparlson  of  Che  3-dlmenslonal-zone- 
folded  Fourier  expansion  calculation  (solid 
line)  for  all  C6X  stage  1  compounds  with  Che 
first  principles  calculation  for  CgLl 
(dotted  line). 3  The  zero  of  energy  is 
ahlfced  to  match  the  degenerate  levels  at 
the  r  point.  The  2s  LI  band  Is  included 
in  the  first  princples  calculation  but  the 
intercalate  band  X  from  the  Fourier 
expansion  calculation  Is  not  included  la 
thn  figure. 


Ce  X  Stage  1 
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Fig.  2  Comparison  of  the  3-dlmensional-zor.e 
folded  Fourier  expansion  calculation  (solid 
line)  for  all  CgX  stage  1  compounds  with  the 
first  principles  calculation  for  CSX  (dotted 
llne).e  The  zero  of  energy  is  shifted  to 
Batch  degenerate  levels  at  the  K  point.  The 
energy  scale  for  CgX  in  this  figure  is  twice 
Chat  of  Ref.  4.  The  4s  and  4p  K  bands  are 
Included  in  the  first  principles  calculation, 
but  the  Intercalate  band  X  from  the  Fourier 
expansion  calculation  is  not  included  in  the 
figure. 
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The  property  of  graphite  Intercalation  coa- 
pounds  that  may  be  of  greatest  commercial  Interest 
la  the  high  electrical  conductivity  resulting  from 
th«  Intercalation  process. ^  Graphite  offers  a 
host  material  with  high  In-plane  mobility  (room 
temperature  mobility  ■  13,000  cm^/V-sec  as  com¬ 
pared  to  Cu  with  35  cm^/V-sec  and  to  Si  with 
1600  cn^/V  -sec), 2, 3  Nevertheless,  graphite  Is  a 
material  with  modest  conductivity  because  of  its 
low  carrier  concentration  (2  x  10~^  carriers/atom 
at  room  temperature)..^  If  we  consider  the  inter¬ 
calate  layer  as  a  source  for  the  Introduction  of 
carriers  Into  'he  graphite  host  (electrons  for 
donor  compounds  and  holes  for  acceptor  compounds), 
then  the  Intercalation  process  provides  a 
mechanism  for  producing  high  conductivity  materials 
by  greatly  Increasing  the  intrinsic  carrier 
concentration  In  a  high  mobility  host  without 
elgnlflcant  loss  of  carrier  mobility.  We  present 
here  a  simple  model  for  the  electrical  conductivity 
which  accounts  for  the  characteristic  features  of 
the  dependence  of  the  in-plane  conductivity  on 
intercalate  concentration.  On  the  basis  of  this 
model  we  can  draw  several  conclusions  about 
characteristics  that  may  lead  to  increased  in-plane 
electrical  conductivity. 

The  simple  phenomenological  model  for  the 
electrical  conductivity  Is  based  on  the  observation 
that  the  total  conductance  per  unit  cell  of  length 
Ic  Is  equal  to  the  sum  of  the  conductances  of  the 
constituent  layers  contained  within  the  unit  cell 
(see  Fig.  1).  In  this  model  we  distinguish 
between  the  graphite  bounding  layers  adjacent  to 
an  Intercalate  layer  and  the  graphite  interior 
layers  that  are  fully  surrounded  by  other  graphite 
layers.  According  to  this  model,  Che  conductance 
equation  Is 

le<W  ■  •‘l<‘’l'®g>  ^  2'o<®gb/‘’g>+<"-2)‘=o<®gl'®g> 

la  which  n  Is  the  stage  of  the  sample,  Co  Is  Che 
graphite  Incetlayer  separation,  and  (di-fcg)  is  the 
separation  of  two  graphite  bounding  layers  between 
which  the  Intercalate  layer  is  sandwiched.  In  a 
conductivity  experiment  we  measure  the  conductivity 
Oa  which  Is  here  conveniently  normalized  to  the  In- 
plane  conductivity  of  pristine  graphite  Og.  Of  the 
various  layers  within  the  unit  cell,  the  conducti¬ 
vity  of  the  graphite  bounding  layers  Ogb  Is 
dominant  because  of  the  high  carrier  density  In 
these  layers  relative  to  the  graphite  interior 
layars^i^  and  because  of  the  ouch  higher  mobility 
of  the  graphite  bounding  layers  relative  to  the 
Incarcalace  layer.  Though  smaller,  the 
COntribuClor.  of  the  graoMte  interior  layers  Cgi 
is  significant,  particularly  for  certain  inter¬ 
calate  species  such  as  the  alkali  metals.  For 
most  Intercalants,  It  is  believed  chat  the 
contribution  from  the  intercalate  layer  la 
negligibly  small  because  the  Fermi  level  generally 
lies  below  the  intercalant  conduction  bands  and 


above  the  Intercalant  valence  levels,^  and  bocause 
ionic  mobilities  tend  to  be  very  low. 

At  low  Intercalate  concentrations,  the 
conductivity  Is  dominated  by  the  fract:  on  of  zne 
unit  cell  occupied  by  the  highly  conduct  i-.g 
graphite  bounding  layers  2co/ (nC(,+d ; ) ,  t.  ^re:  • 
yielding  an  approximately  linear  dependence  o; 
(0,/Og)  on  (1/n)  In  good  agreo:iient  with  much  of  tr.e 
published  experimental  data3>8,9  (see  Fig.  2). 

In  the  low  stage  limit  where  the  semi- 
Insulatlng  intercalate  layer  occupies  a  fignirlccnt 
fraction  of  the  unit  cell,  d^  / (nco^-dj^) ,  izturoticn 
behavior  in  the  (Oa/^g)  vs  Cl/n)  curve  reiul-.i, 
leading  eventually  co  a  decrease  In  (."aZ-r^  Zor 
the  lowest  stage  compounds,  in  agreement  -itr. 
experimental  observations.  This  saturatien  end 
fall-off  effect  at  low  stage  is  more  inpertar.t  for 
Intercalants  with  large  di  values  (c.g.  n!.'03'  and 
less  Important  when  d^  is  small  (e.g.  Li',  also  i- 
sgreemenc  with  experiisencal  observations. 

Contributions  from  the  graphite  interior 
layers  are  significant  when  the  charge  ttir.siTr  to 
these  layers  is  appreciable  and  the  carrier 
mobility  on  the  graphite  interior  layers  us 
significantly  higher  than  on  the  graphite  oouo.dir.g 
layers.  Evidence  for  significant  contribucicr.s 
from  the  graphite  interior  layers  is  the  .-bserva- 
tlon  of  an  enhanced  initial  rise  in  (-a/-:)  vs 
(1/n)  or  alternatively  as  a  shift  in  the  teak  or  the 
onset  of  saturation  of  (Oa/Og)  to  lower  (1/a) 
values. 

Fits  of  this  phenomenological  model  have  beer, 
made  Co  conductivity  data  for  the  donor  ir. tertalur.t 
K  and  the  acceptors  HNOj  and  FeCl3.  Results  of 
this  analysis  show  that  the  conrribution  :ror  the 
graphite  interior  layers  is  more  pronounced  for 
the  K  compounds  than  for  these  acceptor  cu-p:  .r.ds. 
This  conclusion  is  in  good  agreement  with  rtaunetc- 
reflectlon  results  on  the  stage  dependence  o:  the 
Fermi  level  for  the  donor  Rb  and  the  accettots 
FeCl3,  AICI3  and  Brj.lO  Also  of  significunco  is 
the  supporting  evidence  obtained  from  dot s rminot ion 
of  the  dynamic  charge  associated  with  the  orathite 
bounding  and  graphite  interior  layers  trot:  measure¬ 
ment  of  the  strength  of  the  corresponding  in: rare 0 
oscillator  strengths.  Infrared  measuromer. rs  snow 
In  agreement  with  results  on  the  static  c'.arge  that 
the  effective  dynamic  charge  for  the  interior 
layera  la  relatively  larger  for  the  alkali  metal 
donor  compounds  than  for  the  FeCl3  accepter 
compounds. 30 

Thu  phenomenological  conductivity  model 
suggests  that  for  further  enhancement  of  the 
electrical  conductivity  of  graphite  intercalation 
compounds,  it  la  important  to  maximize  the  charge 
density  in  the  graphite  layers  through  increased 
charge  transfer.  Attention  to  Increasing  the 


^Supported  by  the  Air  Force  Office  of  Scientific  Research  (#77-3391). 

^Departaeat  of  Electrical  Englneartng  A  Computer  Science  and  Center  for  Katcrlals  Science  4  Engineering. 
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aaapl*  perfection  for  Incre.iird  carrier  nnblllty  la 
the  graphite  bounding  l.'iver  could  also  be  of 
algntf Icancc.  To  exploit  the  hlch  conductivity  In 
the  graphite  bounding  layers  for  :  iie  low  stage 
coapounds,  the  Intercalate  laver  thickness  dt 
should  be  olniialzed.  The  Interpretation  OL 
conductivity  data  using  this  simple  phenomenolo¬ 
gical  Bodel  thereby  suggests  some  characteristic 
(actors  that  could  be  utilized  in  further 
increasing  the  electrical  conductivity  of  these 
■atarlsls. 
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Ve  present  the  first  detailed  study  of  the  stage  dependence  of  the 
IR-  and  Ranan-aetlve  optic  graphitic  nodes  In  a  graphite  acceptor 
intercalation  compound.  The  general  frequency  upshift  observed  with 
Increasing  FeCl}  concentration  for  all  optic  modes  Is  Interpreted  to 
Indicate  an  In-plsne  compression  within  the  graphitic  layers.  An 
Identlf Icatlon  of  the  IR-activc  modes  with  bounding  and  interior  graphite 
layers  is  siade.  A  lineshape  analysts  of  the  IR  spectra  implies  IR  dipole 
SKiaenta  corresponding  to  'V'lOX  of  the  effective  charge  in  the  graphite 
bounding  layers,  independent  of  stage,  and  'v.JOZ  distributed  among  the 
graphite  Interior  layers  for  stage  n>3  compounds. 


1. INTRODUCTION 

It  is  established  tl-S]  that  the  Raman-active 
lattice  audes  for  graphite  Intercalation  compounds 
can  be  identified  with  (1)  the  intercalate  layer, 
(11)  the  graphite  layers  adjacent  to  the  inter¬ 
calate  layer  (the  graphite  bounding  layers)  or 
(ill)  the  graphite  interior  layers.  Graphitic 
Raman  modes  In  the  vicinity  of  1600  ca~l  have  been 
observed  for  a  variety  of  intercalate  species  as 
doublet  structures  separated  by  approximately 
20  CB“1.  The  lower  frequency  component  (denoted 
by  bSg;)  associated  with  the  interior  gra^^hite 
layersf  and  the  upper  component  (denoted  by  E2g2) 
with  bounding  graphite  layers.  Previous  work 
focussed  on  the  relative  intensities  of  these 
components  as  a  function  of  intercalate  concen¬ 
tration  (1-5]  or  alternatively  as  a  function  of 
stage,  where  the  stage  Index  n  is  the  number  of 
graphite  layers  between  consecutive  intercalate 
layers.  The  present  Raman  study  on  well- 
characterized  graphite-ferric  chloride  samples, 
while  supporting  previous  findings  on  the  inter-  . 
calste  concentration-dependence  of  the  doublet 
separation  and  the  relative  intensity  of  the 
doublet  components,  represents  the  first  systema¬ 
tic  analysis  of  the  graphitic  optic  mode  frequen¬ 
cies  of  an  acceptor  system  over  a  wide  Intercalate 
conecntratlen  range.  The  present  work  also 
represents  the  first  systematic  study  of  the 
Infrared  (IR)  spectra  of  a  graphite  intercalation 
compound.  The  complementary  inforaation  obtained 
from  the  Ran.in  and  Infrared  spectra  has  important 
implications  on  the  interaction  between  intcrc.sl- 
aCe  and  host  layers  in  these  compounds. 

2.FJIPFR1MENTAL  DETAILS 

Vhen  highly  oriented  pyrolytic  graphite  (HOPC) 
Is  employed  as  a  host  material,  staging  conditions 
change  with  respect  to  those  of  single  crystal  or 
flake  graphite  (6-S).  However,  in  agreement  with 


previously  reported  results  for  alkali  metal  [9] 
and  halogen  Intercalants  (10]  it  is  possible  to 
produce  a  series  of  well  staged,  homogeneous 
crystals  over  a  wide  range  of  intercalate  concen¬ 
tration.  We  report  here,  for  the  first  tine, 
results  on  essentially  single  staged  n  •  1,2, 3, 4, 

5,6  and  11  compounds  prepared  from  KOPC  and 
anhydrous  ferric  chloride  (made  in  situ)  using  tlic 
now  conventional  two  zone  system  (11].  Improved 
sample  homogeneity  has  been  achieved  by  u.-^ing 
graphite  of  approximately  20um  initial  thicKness. 

For  the  growth  process  employed,  the  graphite 
temperature  was  maintained  throughout  at  33C*C  and 
the  temperature  (and  hence  vapor  pressure)  of  the 
ferric  chloride  was  varied:  stage  2  (300°C),  stage 
3  (227’C),  stage  4  (215*C),  stage  5  (202*C),  stage 
6  (216*C)  (12)  and  stage  11  (192*C).  For  all  stages 
a  chlorine  atmosphere  of  600  torr  was  employed. 

Using  a  single  crystal  x-ray  diffractometer,  tenven- 
tlonal  0-20  scans  were  taken  to  characterize  the 
samples.  By  employing  a  Si (Li)  detector  and  a 
single  channel  analyzer  to  provide  discrinlnation  of 
incident  photon  energy,  small  concentrations  of 
admixed  stages  were  identified  (13).  X-ray  spectra 
for  stages  n>l,2,3  and  4  are  shown  in  Fig.  1, 
together  with  the  corresponding  intcrcalant  repeat 
distances  I^  determined  using  the  KoKa  x-ray  line. 
The  preparation  of  single  staged  n-1  material  proved 
difficult  because  of  the  tendency  for  inclusion  of 
cither  regions  of  pristine  graphite  or  stage  n>l 
material.  Spectra  for  all  compounds  used  In  this 
study  showed  the  samples  with  stage  n}2  to  be  ciiscn- 
tlally  single-staged.  Only  small  secoudary-onase 
x-ray  diffraction  peaks  (indicated  by  *  in  tig.  i) 
•re  observed.  The  relative  intensities  of  the 
diffraction  peaks  can  be  related  to  the  structure 
factor  for  the  intercalation  compounds  (13|. 

For  the  light  scattering  measurements,  laser 
excitation  was  provided  by  an  argon- Ion  laser 
operating  at  4SS0A  in  the  hack-scac te rt np  geometry. 
Unpolarlzcd  room  temperature  Ranan  spectra  were 
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«».  I.  X-ray  stage  characterization  for  stages 
n  ■  1,  2,  3  and  4  graphite-ferric  chloride 
coapounds.  The  intercalanc  repeat  distance  Ic 
and  stage  indices  are  given  on  the  right.  Reflec¬ 
tions  due  to  adaixed  stages  arc  indicated  by  *. 

The  peaks  labeled  OOtcR  In  (he  stage  n  ••  1  trace 
refer  to  pristine  graphite  reflections.  A  sche- 
•atie  of  the  X-ray  systen  incorporating  a  Si(Li) 
detector  and  a  Ho  X-ray  source  is  shown  in  the 
inset. 


dried  air  and  no  spectral  Intcrforenec  frora  waiter 
vapor  abnorptlon  lines  was  detected  In  the  region 
of  Interest. 

To  check  that  the  reported  x-ray  results  arc 
compatible  with  the  Raman  and  IK  results,  x-ray 
■casurenents  wore  also  made  on  very  thin  samples 
and  on  freshly-cleaved  samples.  Similar  results 
ware  obtained  in  all  cases.  We  have  confidence 
that  our  Raman  and  IR  spectra  are  representative 
of  tho  Indicated  stage. 

3.USIILTS 


Xanaa  spectra  for  graphite  ferric  chloride  are 
shown  In  Fig.  2  for  essentially  single-staged  n  - 
7,3,4,C,11  and  for  our  purest  stage  1  compounds. 
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taksn  from  the  sample  c-faces  and  correspond  to  the 
axeltstlon  of  In-plane  Rooun-setive  modes.  To 
avoid  Intercalate  desorption  effects  associated  with 
laser  heating,  these  spectra  were  taken  at  low 
(<S0bw)  laacr  power  levels. 

Infrared  reflectivity  spectra  were  obtained 
•sing  a  Fourier  transform  spcctrophotoraccer 
operating  with  a  globar  source  and  a  TCS  detector 
vlth  a  KBr  window.  Room  temperature  spectra, 
covering  the  energy  range  dOO  cn~l<u«i000  eo~^  vlth 
O  voeelutlon  of  2  em'i  were  taken  with  the  e-axle 
Opproslm.rcely  parallel  Co  the  beam  path.  The 
opeeirephotomctcc  war  purged  with  an  atmosphere  of 


Fig.  2.  Unpolarizcd  room  temperature  Raman  spectra 
taken  in  the  backscattcrin-  geometry  (E  1  c)  for 
stage  n  •  l,2,3,t,6  and  11  graphite-ferric  chlo¬ 
ride  compounds  and  for  pristine  graphite  (ItOPC). 
Laser  excitation  at  4830.1  and  a  power  level  <  SOmw 
were  used  to  excite  in-plane  Raman-active  F.jgz 
modes  (tec  inset).  The  upper  frequency  component 
(^Jgz)  i*  identified  with  the  bounding  layer  mode 
find  the  lower  component  (E*  jgz)  with  the  interior 
graphitic  mode. 
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'  RaMii  reautca  were  prcvlouely  reported  In  «  pre> 
Itntn.iry  forn  for  moHt  of  these  (14|  and  by 

Caswell  and  Soltii  fur  st.ip.ua  n>l  and  I  IS}. 

'  IneluduJ  for  comp-iriuun  ts  the  spcctrura  for  prts-  ■ 
tine  graplilte  (n""),  showing  a  single  line  at  the 
'it  -  node  frequency.  (See  inset  to  Klg.  2).  Of 
particular  significance  Is  the  large  Inrcrcal.ite 
eonccntratlon-Ucpcndfnt  upshift  oi  the  Ejgj  bound¬ 
ing  layer  code.  The  spectra  for  stages  n"*l  and 
n*2  exhibit  only  a  single  peak,  consistent  with  the 
absence  of  interior  graphite  layers  for  these 
stages.  Saaples  for  n>3  however  show  both  C2gs 
and  C}-.  Bodes,  the  latter  code  associated  wUH 
Interior  graphite  layers.  For  stage  n»4,  where 
the  nunber  of  bounding  and  interior  layers  is 
o^ual,  the  two  peaks  have  approximately  equal 
Intensity.  As  the  intercalate  concentration  de¬ 
creases  and  the  stage  index  Increases,  the  doublet 
separation  reoalns  constant  while  the  relative 
intensities  of  the  Elg,  and  features  exhibit 

n  dependence  on  Intercalate  concentration  (or  1/n) 
in  gualitallve  agreement  with  previous  work  on 
ether  Intercalate  species  (3,4].  Accurate  fits 
to  the  Raman  spectra  have  been  achieved  with  s 
single  Lorentzian  line  for  stages  n-l  and  2  and 
two  Loicntzians  for  stage  n>3.  In  the  lineshape 
analysis,  the  Lorentzlan  lines  were  convolved  with 
the  neasuicd  Instrument  function  and  the  results 
for  the  central  frequency,  llneuldth  and  peak 
intensity  for  the  spectra  in  Fig.  2  are  given  In 
Table  1. 


Infrared  reflectivity  spectra  for  st.ni-r'-  n-T, 
2,6, i  and  11  coapuunds  and  for  gr.'iphltc  .irv  .n 
in  Fig.  3.  IlucauTc  of  i.'ie  liicrca;.cd  f n  r  1.  r 

eoncentraclon  In  the  interc.iSat  luu  conpinuulL;  c<  7\- 
pared  Co  tli.it  of  graphite,  ihc  i. incline  rcdc.t- 
Iv'lty  in  the  vicinity  of  IfcOO  cu”*’  Is  liicto.T..-;'. 
from  ^7GZ  In  graphite  to  •’-fiOZ  for  st.igc  n-ll,  n  1 
to  ‘'"95%  for  stages  n  “  1,2,6  and  6.  Free  c.irri.  r 
effects  are  also  responsible  for  the  major 
in  lineshape  exhibited  by  these  traces.  Throe 
typical  lineshapes  occur,  and  these  are  shoc-n  for 
tiOFC,  stage  n-11  and  stage  n$6.  To  extract  the 
lacclea  node  frequency,  lincwldch  and  oscillator 
atrength,  a  lineshape  analysis  was  performed  and 
results  for  these  parameters  arc  given  In  TaMc  1. 
The  conputed  reflectivity  using  these  parameters 
is  shown  in  Fig.  3  by  the  dotted  curves.  In  this 
analysis,  the  form  for  Che  frce-carrler  and 
electronic  Inccrband  conttibucion.s  Is  largely 
determined  by  the  IR  spectra  observed  over  the 
frequency  range  800<bi<6000  ea~l  (16}.  For  stages 
n«l  and  2,  the  spectra  arc  fitted  with  a  single 
IR-actlve  lattice  mode  at  1186  and  1S31  ca~I, 
respectively.  However,  for  stage  n>'i  two  nodes 
(labeled  and  Ej^)  of  comparable  lincwldth  arc 
required  (sec  Table  1). 

The  dependence  on  reciprocal  stage  (l/n)  of 
the  frequencies  for  the  Raaan-acclve  codes  (Fig.  2} 
and  for  the  IR-active  modes  (Fig.  3}  is  given  in 
Flg;^4.  It  is  noteworthy  that  the  E2g2,  ^2gq*  "lu 
and  mode  frequencies  for  stage  n>3  all  exhibit 


TABU  1  Optic  Mode  Parameters  for  Graphite-Ferric  Chloride  Intercalation  Cocpo’.Tnds 

(Central  frequency  (iJoK  ^nd  strength  (7)  of  Raman-and  IR-ncClve 

nodes)  resulting  from  analysis  of  Che  data  shown  In  Figs.  2  and  3. 


Raman  Infrared 


(a) 

Interior  Laver 

0>) 

Bounding  Layer 

U> 

Interior  Layer 

(b) 

Bounding  Layer 

Stage 

“o 

r  (c) 

r  «) 

R 

u 

0 

r<c) 

y  <d) 

R 

r(c) 

.  (e) 

'iR 

“o 

r(c) 

-  le) 
*IR 

• 

1626.0 

3 

4.0 

1582.5 

8 

0.220 

2 

^612.9 

3 

4.0 

1581.2 

3 

0.110 

3 

1985.0 

13 

3.8 

1607.3 

13 

4.3 

1588.5 

6 

0.025 

1585.8 

4 

0.074 

4 

1583.7 

12 

4.2 

1606.7 

12 

4.2 

1589.4 

6 

0.020 

1585.0 

4 

0.047 

5 

1588.6 

8 

0.015  * 

1584.2 

6 

0.045 

« 

1584.4 

10 

4.2 

1604.9 

9 

2.2 

1588.5 

7 

0.012 

1584.0 

7 

0.037 

9 

1587.3 

4 

0.008 

1583.0 

10 

0.020 

U 

1582.8 

10 

4.4 

1603.8 

6 

2.0 

1587.3 

5 

0.004 

1583.2 

10 

0.020w 

m 

1581.5 

11 

4.8 

1587.5 

2 

0.006^*^ 

(a)  Raman-  and  IR-actlve  graphite  interior  Icyer  modes  are  respectively  denoted  as  E2g2  'lu* 
(k)  Ram^n-  an^  IR-actlve  graphite  bounding  layer  codes  are  respectively  denoted 
••  Ejgj  (R2g2>  a-'l  Eiy(2iu), 

(c)  Central  frequency  uo  and  full-width  half-maximum  lincwldch  T  are  In  cm’*^ 

(d)  Teak  intensities  (arbitrary  units),  accurate  to  ilDZ  except  for  stage  3wtiere‘thc 
uneercalncy  is  n.  t20Z. 

(a)  The  IR  spectra  ate  analyzed  by  c(u)  "  e„  ♦  Fjg  -  «*  -  irw). 

(f)  Value  given  In  Ref.  17  la  'v*  S  tinea  larger. 
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Fi*.  3.  Room  temperature  c-face  infrared  reflectx-  ' 
vity  tpectra  for  the  energy  range  1520$u$l650  cm*‘ 
for  stage  n  «  1,2, 4, 6  and  11  graphite-ferric 
chloride  compounds  and  for  pristine  graphite  (IIOPCl . 
The  fit  CO  experimental  data  using  paraiteters  given 
in  Table  1  is  shovm  by  the  dotted  curves.  These 
data  are  analyzed  in  terms  of  the  graphitic  Eiu 
■odes  (see  inset)  which,  in  the  intercalation 


compounds  are  designated. a 
graphicic  layers  and  as  Eli 
layers:  see  text  and  Table 


a  linear  upshift  in  frequency  with  (1/n)  as  Indl- 
ated  by  the  solid  lines  in  Fig.  4.  Furthernaro, 
all  lines  have  similar  slopes,  resulting  in  doublet 
separations  of  22:2  co”^  for  the  Raman-active  modes 
and  4il  cn~l  for  the  IR-activc  inodes.  Also  shovm 
In  Fig.  4  (see  inset)  is  Che  stage  dependence  of 
<Fa>,  the  In-planc  dipole  moment  for  the  modes 
resulting  from  the  lineshape  analysis  of  the  IR 
data.  In  this  inset,  Che  value  of  <Pa>  is 
notmalised  (i.e.  1  as  n  *  *). 

4.D1SCUSSI0.Y 

The  main  eonttlbutions  of  this  work  are:  (1) 
identification  of  the  observed  iK-aetive  nodes 
vith  interior  and  bounding  graphite  l.iyers,  (2) 
determln.K ion  of  the  rliarge  assorlated  with  the 
dipole  twcK'nt  in  the  bounding  ,ind  interior  layers. 


Eiu  for  anterior 
I  for  bounding  graphitic 

1. 


and  (3)  observation  of  a  similar  stage-dependent 
upshift  in  frequency  for  all  Raaan-actlvc  and  IR- 
active  gr.iphitic  modes  for  stage  n>3  compounds. 

The  IR  mode  labeled  E^u  in  Fig.  4  is  identified 
with  vibrations  in  the  graphite  interior  layers 
for  the  following' reasons.  Firstly,  the  frequency 
of  this  mode  extrapolates  to  that  for  the  Exu  "odc 
for  graphite  as  the  intercalate  concentration  goes 
to  tcro  (ri'e°i).  Secondly,  this  mode  is  not  found 
for  stage  n>l  or  2  compounds  (where  there  arc  no 
interior  graphite  layers),  but  is  observed  only 
for  stage  nA3  compovmds,  FuTthemore,  analysis 
of  the  IR  lincshapcs  shows  that  the  dipole  moment 
<Pa>  (see  inset  to  Fig.  4)  for  the  E{„  node  in  tlic 
Interc.alatlon  compounds  decreases  with  tncrcasing 
stage  index  and  excr.tpolatcs  to  the  v.tliiu  for 
grophltc.  Consistent  with  this  Idcncific.ation  for 
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FIr.  4.  The  Inverse  stage  dependence  of  the  node 
frequencies  for  Raman  E’g2  and  Ejgj  (squares)  and 
IK  E,u  and  E^u  <c  ircles)  modes  in  graphite-ferric 
chloride  compounds.  The  solid  lines  represent  a. 
least  square  fit  to  the  cxperincntal  points  for 
li>3  and  indicate  a  similar  strain  dependence  for 
all  observed  modes.  The  and  Eiu  mode  frequencies 

for  pristine  graphite  are  given  as  (l/n)->0.  The 
inset  shows  the  (1/n)  dependence^of  the  dipole 
swincnt  for  the  infrared  E^u  and  Eiu  modes  normalized 
to  that  for  the  Eju  mode  in  pristine  graphite.  The 
KesHn  and  IR  frequencies  for  the  lowest, stages 
B«l  and  n»2  compounds  labeled  £.252  and  Eju  depart 
from  the  behavior  observed  for  the  ni3  compounds. 


the  Eiu  mode  Is  the  observation  that  the  magnitude 
of  <Pa*  for  this  IR  mode  incrc<ascs  by  a  factor  of 
four  from  its  small  value  in  graphite  (17). 

A  second, IR-actlvc  mode,  absent  in  graphite 
and  labeled  Ej^  in  Fig.  4,  is  identified  witii  the 
bounding  layer.  The  dipole  moment  for  this  m.udc 
has  a  much  weaker  stage  liependence  and  tends  to 
ba  large  compared  with  <Pa>  for  the  eJu  mode  for 
high  stage  compounds  (see  inset  to  Fig.  4).  ^hc 
approximate  equality  of  <Pa>  for  the  eJu  and  Eiu 
modes  at  stage  3  is  also  consistent  with  this 
identification.  The  single  IR  mode  for  the 
lewect  stage  n(2  compounds  is  clearly  due  to 
hounding  layer  vlhractons.  The  frequencies  of  the 
hounding  l.sycr  nodes  arc  anomalously  low  for  stages 
n“l  and  n“2,  while  the  corresponding  frequencies 
of  the  Rasitin-Act ivc  nodes  are  anomalously  lilgh. 

Por  stage  2,  the  bounding  layer  lies  between 
another  bounding  layer  and  an  Intercalate  layer. 


while  for  stage  1  the  bounding  layer  lies  between 
two  intercalate  layers.  In  both  cases  therefore, 
the  intcrplanar  coupling  Is  different  from  the 
situation  for  stage  n>3  where  the  bounding  l.aycr 
is  adjacent  to  an  interior  graphite  and  an  inter¬ 
calate  layer.  Analysis  of  the  IR-modc  intensity 
provides  values  for  the  effective  charge  e* 
associated  with  the  dloolc  moment  [17).  Since  the 
and  occur  at  different  frequencies,  the 
relative  charge  associated  with  the  graphite 
interior  and  bounding  layers  (18,19)  can  be 
deCCTmlned  Independently.  The  data  In  Table  I 
were  analyzed  using  the  Szigett  effective  charge 
model  (20),  assuming  that  per  unit  cell  there  Is 
one  osclll.ator  per  hounding  layer  and  (n-2) 
oscillators  dlstribuced  among  the  interior  lavers. 
With  these  .istumpt  lens  tlie  data  yield  rea  -on.iele 
values  for  the  •effective  charge  c*  per  unit  c.  ll; 
a*  ^  0.73c  for  stages  n*l  and  2  and  e*  s*  0.93e  for 
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n>3.  la  adUltloa,  aaolyKls  shows  Lhac  on 
•ach  boundlnc  r.<>>!'hlcc  layer  ts  ap;iroxln;icclv  !:L.'ip,e 
Independent.  Kor  sta?.c  n>l,  the  cf<i-ctl«o  charr.e 
distributed  anonr;  the  Interior  rinphltc  layers  Is 
nsscntlally  constant  at  ICiZ,  with  the  re.-r-ilnlng 
70Z  residing  on  the  bounding  layers.  Assunlng 
the  sane  charge  In  each  granliite  Interior  layer 
then  the  charge  per  Interior  layer  Is  stage 
dependent.  .. 

The  upshift  of  the  Eiu  and  Elu  nodes  with 
Increasing  Intercalate  concentration  for  stage 
n^  Is  in  good  agrecaent  with  the  upshift  observed 
for  the  Ejg2  ^2g2  Raraan-actlvc  nodes.  Ue 
attribute  the  stiffening  of  these  lattice  nodes  to 
n  strain  which  decreases  the  in-plane  graphite 
lattlec  constant  [21).  The  eonstancy  of  the  line- 
width  for  the  F.222  interior  layer  code  as  a 
function  of  intercalate  concentration  and  an 
Identical  frequency  upshift  for  Loth  Interior  and 
bounding  Ranan  anJ  IR-activc  nodes  suggests  that 
on  Identical  strain  exists  within  both  the  bound¬ 
ing  and  Interior  graphitic  layers  for  a  given 
coopound.  On  the  other  hand,  the  oagnltude  of  the 
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strain  Incre.asea  as  the  number  of  Interior  graplilto 
layers  i)eer>'ases,  cen.slstcnt  with  the  idea  rh.u  the 
stress  Introduced  by  Intercalation  most  be  stmed 
by  fewer  layers  as  n  decreases  (7?1.  Pic  coi  ii;:- 
dence  In  the  Ej„  and  Ei..,  node  f  rcnucnc  les ,  the 
relatively  snail  lincwidtii  ot  the  K.anan-acc ive 
bounding  layer  node,  and  the  anonalous  .hanin  and  IK 
awde  frequencies  for  stage  n*l  and  2  conpounds 
present  Important  characteristics  which  must  be 
explained  by  a  theory  for  the  lattice  mode  spectra 
of  graphite  Intercalation  coapounds. 
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ABSTRACT 

A  theoretical  model  for  the  electronic  structure  of  the  graphite 
intercalation  compounds  based  on  the  k  axis  zone  folding  of  the 
pristine  graphite  ir-bands  is  applied  to  explain  the  experimental 
results  obtained  in  magneto-reflection  and  in  Shubnikov-de  Haas 
experiments  on  donor  compounds.  The  detailed  fit  to  the  observed 
Shubnikov-de  Haas  frequencies  for  stage  5  graphite-potassium  supports 
the  model  and  suggests  that  the  dominant  interactions  that 
determine  the  Fermi  surfaces  are  closely  related  to  those  involving 
the  ir-bands  for  pristine  graphite. 
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I.  Introduction 

The  recent  development  of  theoretical  models  for  the  electronic 
structure  of  graphite  Intercalation  compounds  gives  promise  for  new 
advances  in  our  understanding  of  this  Important  class  of  compounds. 

These  models  include  first  principles  calculations  for  stage  1  alkali 
metal  compounds  [1,2],  modified  two  dimensional  graphite  layer  models  [3], 
and  models  based  on  the  three  dimensional  graphitic  Slonczewski-Weiss-McClure 
(SUMcC)  model  using  the  band  parameters  for  pristine  graphite  [4]. 

A  common  feature  of  all  of  these  approaches  is  that  the  ir-bands 
derived  from  pristine  graphite  are  the  dominant  electronic  states  in  the 
Intercalation  compound  within  an  eV  of  the  Fermi  level.  This  paper  focusses 
on  the  models  which  are  directly  related  to  the  three-dimensional  pristine 
graphite  model  Hamiltonian  of  SUMcC  and  demonstrates  that  a  phenomenological 
energy  band  structure  can  be  developed  which  is  consistent  with  a  number  of 
experimental  magneto-optical  and  magneto-transport  measurements  on  these 
compounds.  Since  a  large  number  of  Intercalation  compounds  have  been 
prepared,  much  work  remains  to  be  done  to  fully  establish  the  parameters 
of  the  present  model  or  Indeed  whether  some  alternate  model  gives  a  more 
convergent  treatment  for  a  particular  class  of  graphite  intercalation  compounds. 

The  c-axls  zone  folded  model  considered  here  [4]  is  a  simple  extension 
of  the  dilute  limit  model  proposed  earlier  and  used  by  Suematsu  et  al.  [5] 
to  interpret  de  Haas-van  Alphen  (DHVA)  data  and  by  Mendez  et  al.  [6]  to 
Interpret  magnetoreflection  data.  Since  these  two  experiments  played  a 
dominant  role  in  the  evaluation  of  the  SWMcC  parameters  in  pristine  graphite, 
they  will  be  the  experiments  which  are  the  focus  of  this  paper  for  the 
intercalation  compounds.  This  paper  will  show  that  without  introduction  of 
additional  parameters  (i.e.  the  intercalate  layer  is  approximated  as  an 
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empty  layer  between  two  graphite  bounding  layers)  a  model  Hamiltonian  Is 
developed  which  accounts  for  the  main  features  of  the  Shubnikov-de  Haas 
(SdH)  and  magnetoreflection  spectra  obtained  for  alkali  metal  donor  compounds. 

II .  Band  Model 

To  obtain  electronic  dispersion  relations  E(^)  for  graphite 
intercalation  compounds,  a  Hamiltonian  is  developed  based  on  the  SVMcC 
three-dimensional  model  for  the  graphite  7r-bands.  As  a  first  step,  the 
c-axls  superlattice  periodicity  is  included  and  an  "empty  intercalate 
lattice"  for  the  intercalation  compounds  is  Introduced,  depending  only 
on  the  graphite  band  parameters  that  are  already  known.  Intercalant- 
speclflc  Interactions  between  the  Intercalant  and  the  graphite  bounding 
layer  can  then  be  introduced  to  obtain  the  final  dispersion  relations. 

The  main  features  of  this  model,  presented  elsewhere  [4],  are  summarzled 
below. 

The  4x4  SWMcC  Hamiltonian  for  the  graphite  Tr-bands  H  (k  )  is  zone 

o  '  z 

folded  along  k  to  account  for  the  staging  periodicity  as 
z 


Hr  ,  -  ) 

folded  s 


O  8 


(1) 


in  which  k  is  the  wave  vector  near  the  HK  axis,  each  of  the  blocks  is  the  4x4 

9 

SWMcC  Hamiltonian  and  1  *  n+l  for  even  stage  or  (n-M)/2  for  odd  stage  compounds, 
where  n  is  the  stage  index. 

r  In  order  to  consider  explicitly  the  effect  of  intercalation,  the 
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Hamiltonian  given  by  Eq.  (1)  must  be  transformed  to  a  layer  representation 
which  is  written  explicitly  by 


H.  (k  ) 
layer  s' 


U  H 


folded 


rl 


(2) 


where  the  unitary  transformation  is  derived  from  the  basis  functions,  yielding 
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The  Intercalation  process  is  modelled  by  substitution  of  an  Intercalate  layer 
for  one  or  more  of  the  graphite  layers  in  the  layer  Hamiltonian.  The 
"empty  Intercalate  lattice"  model  for  a  graphite  intercalation  compound  is 
then  obtained  by  treating  the  intercalate  layer  as  an  empty  layer,  neglecting 
interactions  within  the  layer  and  between  this  layer  and  the  other  (21-1) 
layers.  The  "empty  intercalate  lattice"  model  thus  requires  no  band  parameters 
other  than  those  for  graphite  which  are  already  well  established  [7].  To 
obtain  the  full  model  for  the  dispersion  relations.  Interactions  between  the 
Intercalate  layer  and  the  graphite  bounding  layers  are  introduced  and  the 
pertinent  band  parameters  evaluated  by  a  fit  to  experimental  data.  In  the 
present  work  the  two  experiments  that  are  emphasized  are  the  Shubnlkov-de 
Haas  and  magnetoreflection  experiments  which  are  respectively  sensitive  to 
the  Fermi  surface  and  the  energy  band  structure  near  the  Fermi  level. 

In  Fig.  1  results  for  E(k)  for  the  "empty  intercalate ‘lattice" 
model  for  stage  1,2  and  3  graphite  Intercalation  compounds  are  given.  These 
''empty  intercalate  lattice"  bands  are  derived  for  an  in-plane  2x2  superlattice 
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and  for  no  intercalace-graphice  bounding  layer  InCeracclon.  Only  Che 
graphite  ir-bands  are  shovm  and  these  bands  are  also  plotted  on  an  expanded 
scale  along  and  near  the  H-K  axis  of  the  folded  hexagonal  Brillouin  zone. 

For  this  empty  Intercalate  lattice  model,  the  TT-bands  near  the  H-K  axis 
are  Independent  of  in-plane  folding  and  in  the  limit  of  stage  these 
dispersion  curves  reduce  exactly  to  the  SWMcC  ir-bands.  Comparison  of  the 
stage  1  results  [8]  with  the  band  structure  of  Inoshlca  et  al.  [9]  gives 
good  agreement  for  the  graphitic  bands  between  the  two  methods .  The  H  and 
K  points  in  the  Brillouin  zone  are  high  symmetry  points  and  show  extremal  cross 
sectional  areas  which  are  expected  to  give  rise  to  de  Haas-van  Alphen 
oscillations.  Since  the  H  and  K  points  are  in  general  non-degenerate,  a 
three-dimensional  calculation  is  needed  to  obtain  different  DHVA  frequencies 
for  K  and  H  point  extremal  cross  sectional  areas. 

III.  Magnetoref lection  Results 

The  magnetoreflection  experiment  is  important  in  providing  detailed 
information  on  the  form  of  E(^)  for  the  nearly  degenerate  bands  near  the 
K-point  in  the  Brillouin  zone  (see  Fig.  1).  In  the  magnetoreflection 
experiment,  resonant  behavior  is  observed  in  the  magnetic  field  dependence 
of  Che  reflectivity  and  these  resonances  are  identified  with  Landau  level 
transitions  associated  with  these  nearly  degenerate  bands  from  occupied 
valence  magnetic  subbands  to  unoccupied  conduction  subbands  above  the  Fermi 
level.  To  relate  these  results  to  band  models,  the  experiments  are  carried 
out  on  well-staged  samples  and  changes  in  the  magnetoreflection  spectra 
are  studied  as  a  function  of  stage  and  intercalant.  Typical  magnetoreflection 
results  [6]  are  given  in  Fig.  2  for  an  acceptor  compound  (stage  7  graphite- 
FeClj)  and  for  a  donor  compound  (stage  6  graphite-Rb) .  Resonances  identified 
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with  R-point  Landau  level  transitions  are  indicated  by  arrows  and  the 
quantum  numbers  for  the  initial  valence  and  final  conduction  magnetic 
aubbands  are  given  in  the  figure.  Experimental  results  on  a  variety  of 
acceptor  and  donor  compounds  Indicate  that  the  form  for  the  dispersion 
relations  E(k)  for  bands  close  to  the  K-point  degeneracy  remains  similar 
to  that  in  graphite,  consistent  with  the  model  for  E(k)  presented  in  Section  II. 
On  the  other  hand,  small  but  measurable  changes  in  the  graphite  band 
parameters  are  obtained  and  evaluated  quantitatively.  Results  for  several 

representative  compounds  are  given  in  Table  1  in  terms  of  valence  and 

*  * 

conduction  band  effective  masses  m'  and  m  and  the  ratio  of  SWMcC  band 

V  c 

2 

parameters  (y  /Y.)  which  sensitively  determine  the  Landau  level  separation, 
o  i 

Furthermore,  analysis  of  magnetoreflection  structure  which  in  the  case  of 
pristine  graphite  is  associated  with  H-polnt  Landau  level  transitions 
indicates  that  y^,  the  nearest  neighbor  overlap  Integral  is  changed  by  less 
than  5Z  for  all  the  Intercalate  species  and  stages  that  have  been  studied  [6]. 
Since  the  S'.VtIcC  band  parameters  are  not  strongly  perturbed  by  Intercalation, 
the  fit  to  the  Fermi  surface  measurements,  discussed  in  the  next  section, 
is  made  as  a  first  approximation  assuming  the  full  set  of  SUMcC  parameters 
appropriate  to  pristine  graphite. 

IV.  Shubnikov-de  Haas  Results 

Study  of  the  Shubnlkov-de  Haas  effect  in  graphite  Intercalation 
compounds  provides  information  on  the  Fermi  surfaces  associated  with  each 
of  the  partially  occupied  bands.  To  obtain  quantitative  information  on 
the  stage  dependence  of  these  Fermi  surfaces  it  is  important  to  work  with 
well-staged  and  characterized  samples.  Results  obtained  for  the  alkali 
■etel  compounds  with  the  Intercalants  K  and  Rb  show  stage-dependent  Fermi 
•urfecas,  so  that  staging  fidelity  was  essential  for  obtaining  reproducible 


experimental  results  for  different  samples  of  the  same  stage  index. 

The  samples  were  prepared  using  the  conventional  two-zone  technique 
with  HOPG  (highly  oriented  pyrolytic  graphite)  as  the  host  material. 

The  sample  dimensions  were  typically  (5x6x0. 5  mm).  Because  of  the  instability 
of  alkali-metal  samples  in  the  presence  of  air  and  moisture,  the  samples 
were  encapsulated  in  ampoules  and  sample  handling  was  done  in  an  Argon- 
filled  dry  box  (M  ppm  oxygen  content) .  The  stage  of  the  samples  was 
determined  using  (001)  x-ray  diffraction  profiles  both  before  and  after 
the  SdH  experiments,  confirming  that  the  samples  were  single-staged  and  that 
no  desorption  had  occurred  during  the  measurements. 

The  four-point  method  was  used  to  study  the  a-plane  transverse  magneto- 
resistance  in  the  temperature  range  1.4  <  T  <  4.2K  and  in  magnetic  fields 
up  to  IS  Tesla.  The  leads  were  attached  to  the  sample  using  conducting 
epoxy.  The  sample  was  then  Inserted  in  a  helium-filled  ampoule  and  stycast 
was  used  to  seal  the  ampoule.  The  angular  dependence  of  the  SdH  oscillations 
could  thus  be  measured  by  rotating  the  sample  around  the  direction  of  the 
current  I  such  that  I  H  for  all  angles,  and  hence  transverse  magneto- 
resistance  was  always  measured.  Data  acquisition  was  by  computer,  and  the 
data  were  manipulated  to  obtain  a  Fourier  power  spectrum  of  resistance 
vs.  1/H,  thereby  yielding  the  frequencies  of  SdH  oscillations.  These 
frequencies  are  related  to  the  extremal  cross  sections  of  the  Fermi  surface. 
The  Fourier  power  spectra  for  the  SdH  frequencies  for  Rb  stages  2,8  and  K 
stages  2,5,8  are  presented  in  Figs.  3  and  4,  respectively. 

During  the  course  of  these  studies,  a  great  deal  of  care  was  given  to 
ensure  the  fidelity  and  reproducibility  of  the  data.  For  this  reason,  the 
experiment  was  performed  on  two  potassium  stage  5  samples  and  on  one  of 
these,  the  measurements  were  repeated  several  weeks  later  after  attaching 


new  leads.  While  x-ray  profiles  after  each  step  of  the  experiment  confirmed 
that  the  samples  retained  their  single  stage  identity,  the  magnetoresistance 
oscillations  and  the  relevant  Fourier  power  spectra  revealed  identical 
traces  for  these  stage  5  samples.  Measurements  done  on  potassium  samples 
with  different  stages  or  on  rubidium  samples,  however,  showed  distinctly 
different  SdH  frequencies.  Hence  our  conclusion  is  chat  there  is  a  unique 
set  of  SdH  frequencies  associated  with  each  stage  and  donor  intercalanC  (K,Rb). 
Suematsu  et  al.  [5]  have  also  shown  stage  dependent  de  Haas-van  Alphen  frequencies 
In  stages  3  and  4  potassium  compounds. 

The  model  for  E(k)  described  above  is  now  applied  to  the  interpretation 
of  the  observed  SdH  frequencies.  Specific  application  is  here  made  for  stage 
S  graphiCe-K,  and  in  this  connection  the  energy  band  model  for  n»S  is  shown 
in  Fig.  5,  for  the  five  valence  and  five  conduction  ir-bands.  In  order  to 
calculate  the  SdH  frequencies  Che  Fermi  level  must  be  determined.  The 
empirical  position  of  Ep  shown  in  Fig.  5  yields  four  partially  occupied 
conduction  bands  with  the  Fermi  surface  parameters  listed  in  Table  2, 

Including  calculated  and  observed  SdH  frequencies,  the  cyclotron  mass  at  Ep,  the 
trigonal  warping  anisotropy  listed  as  and  the  electron  density  for 

each  carrier  pocket.  The  values  for  2  suggest  that  trigonal 

warping  Is  Important  for  the  larger  cross-sectional  areas  with  heavier  masses, 
whereas  the  smaller  light  mass  cross-sections  are  circular.  The  generally 
good  agreement  of  Che  observed  SdH  frequencies  with  the  empty  intercalate 
leyer  swdel  of  Fig,  5  suggests  that  the  effect  of  intercalant-graphlte 
bounding  layer  interactions  can  be  treated  as  a  perturbation  and  evaluated 
by  fitting  the  model  quantitatively  to  the  observed  SdH  frequencies.  It  should 
be  noted  that  the  volumes  of  the  carrier  pockets  for  the  four  occupied  bands  in 
Fig.  S  correspond  to  a  charge  transfer  of  M).3  electron  per  Intercalant  into 
these  carrier  pockets. 
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V.  Concluding  Remarks 

The  nathematlcal  technique  used  here  to  obtain  the  electronic  dispersion 
relations  can  also  be  applied  to  obtain  phonon  dispersion  relations  for 
graphite  Intercalation  coopounds  [10].  A  aajor  strength  of  this  general 
approach  is  the  proper  treatment  of  symmetry  In  the  Fourier  expansion  of 
the  dispersion  relations  for  the  Intercalation  compounds. 

A  three-dimensional  Fourier  expansion  of  the  electronic  energy  bands 
for  intercalated  graphite  based  on  a  k^-axis  folded  SUMcC  structure  has  been 
developed.  An  "empty  intercalate  layer"  limit  of  this  model  has  been  applied 
to  Interpret  the  observed  SdH  frequencies  in  stage  5  graphite-potassium. 

Using  only  the  SNMcC  parameters  for  pristine  graphite,  and  no  other 
parameters,  the  model  Is  able  to  account  for  most  of  the  observed  SdH 
frequencies.  The  good  agreement  with  the  magnetoreflection  results  Is  built 
Into  the  model  by  using  the  basic  SUMcC  Hamiltonian  In  which  the  band  parameters 
were  originally  evaluated  to  fit  the  observed  magnetoreflection  spectra  for 
graphite.  Further  tests  of  the  k  -axis  folded  SUMcC  model  are  now  in  progress 
by  applying  the  model  to  the  SdH  spectra  for  other  stages  in  the  graphite- 
potassium  system  and  for  compounds  with  the  intercalant  Rb.  The  model  Is 
also  being  tested  by  an  application  to  the  Interpretation  of  measured 
optical  reflectivity  spectra. 
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Fig.  1 


Fig.  2 


Fig.  3 


Fig.  4 


Fig.  5 


FIGURE  CAPTIONS 

Electronic  energy  levels  derived  by  k2‘'axls  zone-folding 
of  Che  3  dimensional  Fourier  expansion  of  Che  pristine 
graphite  ir-band  Hamiltonian  for  a  primitive  (2x2)  superlactice.  An 
empty  intercalate  layer  is  assumed  and  the  expansion  parameters 
are  based  on  the  SVMcC  parameters  along  the  HK  axis.  On  the 
right  an  expanded  scale  is  used  to  plot  the  levels  on  and  near 
the  HK  axis. 

MagneCoreflectlon  spectrum  using  (-)■)  circular  polarization  for 
an  acceptor  compound  (FeCl^  stage  7)  at  a  photon  energy 
dSiu  *  0.295  eV  and  for  a  donor  compound  (Rb  stage  6)  at 
dim  *  0.335  eV.  For  comparison,  traces  for  graphite  are 
shown  at  comparable  photon  energies.  The  resonances  are 
specified  by  the  quantum  numbers  for  the  Initial  and  final 
states.  (From  the  work  of  Mendez  et  al.  Ref.  6). 

Shubnikov-de  Haas  Fourier  Transform  Power  Spectra  for  stage  2 
and  8  graphite-Rb.  These  power  spectra  were  obtained  by  a 
Fourier  transform  of  an  experimental  resistivity  vs  1/H  trace 
for  magnetic  fields  H  <  15  Tesla.  The  peaks  in  the  power  spectra 
correspond  to  SdH  frequencies,  which  are  given  in  Tesla  and  the 
same  scale  is  used  for  each  stage. 

Shubnikov-de  Haas  Fourier  Transform  Power  Spectra  for  stage  2, 

5  and  8  graphite-K.  (See  Fig.  3). 

Electronic  band  structure  in  the  vicinity  of  the  HK  axis  near  the 
Fermi  level  for  a  stage  5  "empty  Intercalate  layer"  model  for  a 
graphite  intercalation  compound.  The  Fermi  level  is  determined  to 
fit  the  observed  SdH  frequencies.  For  the  indicated  Fermi  level  four 
of  the  conduction  bands  are  partially  occupied  and  give  rise  to 
Fermi  surfaces  and  carrier  pockets. 


Table  2  Feral  Surface  Parameters  Associated  with  Stage  5  Graphlte^Potasslua 
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PHENOMENOLOGICAL  MODEL  FOR  THE  ELECTRONIC  STRUCTURE 


OF  GRAPHITE  INTERCALATION  COMPOUNDS 
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ABSTRACT 


The  electron  energy  bands  for  graphite  intercalation 
compounds  are  calculated  uslnlg  a  model  based  on  the  well- 
established  electronic  structure  of  pristine  graphite.  The  ir-band 
Hamiltonian  for  graphite  is  folded  in  k-space  appropriate  to  the 
Intercalant  superlattice  which  includes  k^-axls  and  in-plane  zone 
folding.  The  folded  k-space  Hamiltonian  is  transformed  into  a 
layer  representation  and  modifications  are  incorporated  into  the 
Hamiltonian  to  properly  account  for  the  perturbations  due  to  the 
intercalate  layer.  With  the  inclusion  of  an  interaction  between 
the  intercalate  layer  and  the  graphite  bounding  layer,  electronic 
structures  appropriate  to  both  donor  and  acceptor  compounds  can 
be  deduced.  The  model  is  applicable  to  any  intercalant  and  stage, 
and  predicts  stage  dependent  Shubnikov-de  Haas  and  optical  spectra 


1.  Introduction 


The  intercalate  symmetry  imposes  a  superlattice  periodicity  on  the 
crystal  structure  of  graphite  intercalation  compounds  which  is  reflected 
In  the  electronic  structure.  This  Intercalate  symmetry  can  be  conveniently 
modeled  by  Brillouln  zone  folding  of  the  electronic  dispersion  relations 
E(k)  of  pristine  graphite  [1].  The  staging  periodicity  is  the  dominant 
symmetry  of  Intercalated  graphite,  and  can  be  found  in  most  donor  and  acceptor 
compounds  that  have  been  prepared  to  date  [2].  The  staging  periodicity  is 
treated  mathematically  in  terms  of  a  k^-axis  zone  folding  of  E(k).  In 
addition,  in-plane  ordering  appears  to  be  well  established  in  some  compounds, 
such  as  the  stage  1  alkali  metal  compounds  [3].  The  symmetry  imposed  by  the 
In-plane  Intercalate  ordering  for  systems  where  the  intercalate  structure  is 
commensurate  with  the  graphite  layers  is  treated  by  an  In-plane  folding  of 
the  dispersion  relations.  Both  types  of  zone  folding  are  considered  in 
this  treatment  of  the  electronic  dispersion  relations  for  intercalated 
graphite.  By  emphasizing  the  basic  symmetry  of  the  Intercalation  compounds, 
a  model  is  developed  which  is  applicable  to  any  intercalant  and  stage. 

The  zone  folding  technique  presented  here  is  based  on  the  known  electronic 
structure  of  pristine  graphite  [4,5].  The  intercalate  layer  is  initially 
modeled  as  a  periodic  empty  layer  [6]  in  otherwise  perfectly  ordered  graphite. 
This  empty  intercalate  layer  model  simultaneously  explains  Fermi  surface 
data  [7],  which  show  a  pronounced  stage  dependence  of  the  Fermi  surfaces 
and  major  departures  from  the  graphite  Fermi  surfaces,  and  magnetoref  lection 
data  [8]  which  provide  .strong  evidence  that  the  enetgy  bands  near  the  Fermi 
level  remain  highly  graphitic  upon  Intercalation.  A  k  -dependent  charge 
distribution  results  from  the  incomplete  occupation  of  zone  edge  valence 
or  conduction  bands  and  from  the  Intercalant-graphite  bounding  layer 
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Interactions  which  are  introduced  to  obtain  a  quantitative  fit  to 
experimental  data  . 

In  this  paper,  explicit  results  for  the  electronic  dispersion 
relations  are  presented  for  various  In-plane  superlattlces  for  stage  1 
compounds,  and  for  some  higher  stage  compounds.  Calculated  dispersion 
relations  are  also  presented  as  a  function  of  stage,  and  are  related  to 
experimental  results  relevant  to  the  electronic  structure. 


II.  Brillouin  Zone  Folding  Model 

The  introduction  of  the  intercalant  introduces  a  c-axis  periodicity 
associated  with  the  staging  phenomenon  and  in  some  cases  an  additional 
in-plane  periodicity  when  the  intercalant  is  commensurate  with  the  graphite 
In-plane  structure.  The  effect  of  the  superlattice  periodicity  on  the  electronic 
structure  is  treated  by  zone  folding  of  the  graphite  Brillouin  zone  on  to  the 
smaller  zone  for  the  superlattice  structure.  The  Interaction  between  the 
Intercalate  and  graphite  is  then  treated  in  perturbation  theory. 

In  cases  where  an  in-plane  superlattice  is  present,  the  zone  folding 
is  represented  mathematically  by  the  in-plane  zone-folded  Hamiltonian 


H  (k+K.)  . 
o  -  -X 


H  (k+K  ) 
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(1) 


in  which  H^(k)  is  the  4x4  Hamiltonian  for  the  graphite  n-bands  and  the 

reciprocal  lattice  vectors  K, ...K  are  given  in  Table  I  for  the  simpler 

.X  _m 
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In-planc  structures.  For  a  superlattice  unit  cell  containing  p  graphite 
unit  cells  there  are  p  reciprocal  lattice  vectors  and  Hjy^(k)  is  a  (4px4p) 
■atrix.  If  there  is  no  in-plane  periodicity,  then  Hj^^(k)  is  simply  the 
(4x4)  Slonczewskl-Welss-McClure  or  Fourier  expansion  Hamiltonian  for 
pristine  graphite. 
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The  staging  periodicity  for  the  Intercalation  compound  Is  modeled  using 

a  k.  -axis  zone  folding  technique.  The  Hamiltonian  is  written  in  terms  of 
z 

£  pairs  of  graphite  planes  where  the  basic  repeat  distance  is 

determined  by  the  number  of  planes  in  the  c-axis  repeat  distance  of  the 
Intercalation  compound,  considering  for  the  moment  the  intercalate  layer 

to  be  indistinguishable  from  graphite.  The  k  -zone-folded  Hamiltonian 

z 

H^^(k)  is  then  written  as 

H|^(k)  0  ... 

H  M  -  “  ’’if  *  IT  (2) 

ZZ  ~  O 

•  • 

« 

"If  <!;  *  —  - 
o 

where  each  of  the  blocks  is  the  (4px4p)  in-plane  folded  graphite  Hamiltonian 
given  by  Eq.  (1)  and  I  =  n+1  for  even  stage  or  (n+l)/2  for  odd  stage  compounds, 

O 

where  n  is  the  stage  index  and  c^  =  3,35A  [1], 

The  effect  of  intercalation  is  easily  implemented  in  the  Hamiltonian 
given  by  Eq.  (3)  after  making  a  transformation  to  a  layer  representation  [1]^ 
which  is  written  as 

where  the  unitary  transformation  transforms  the  basis  functions  of  Eq.  (2) 
to  basis  functions  for  a  layer  representation.  The  resulting  layer 
Hamiltonian  is  then  written  as 


in  which  the  blocks  2px2p  matrices  having 

layer  subscripts. 

The  unitary  trans forma cion  of  H^^(k)  does  not  change  the  eigenvalues 
of  the  Hamiltonian.  Thus  the  energy  bands  from  the  folded  graphite 
Hamiltonian  in  the  layer  representation  of  Eq.  (A)  are  formally  identical 
to  the  graphite  energy  bands.  However,  in  the  layer  representation,  the 
effect  of  intercalation  is  easily  included  by  substitution  of  matrix  blocks 
for  the  intercalate  species  in  place  of  the  appropriate  graphite  layer. 

For  the  odd  stage  n  compounds,  the  resulting  Hamiltonian  is 


where  \  is  written  for  (n-l)/2,  and  X  denotes  an  intercalate  layer.  For 


the  even  stage  compounds  [1]  both  an  A  layer  and  a  B  layer  must  be  replaced 
by  an  intercalate  block.  Explicit  results  for  the  electronic  energy 
dispersion  relations  are  now  presented  for  several  cases  of  Interest. 

III.  Results  and  Discussion 

Dispersion  relations  for  graphite  and  for  various  in-plane  zone  foldings, 
some  of  which  are  appropriate  to  stage  1  compounds  in  their  low  temperature 
phases, are  shovm  in  Fig.  1.  These  TT-band  dispersion  relations  are  calculated 
from  Eq.  (5)  for  the  case  where  Che  graphite  layer  is  replaced  by  an  empty 
intercalate  layer  (the  "empty  Intercalant  layer  lattice") .  We  note  that  in 
constructing  the  dispersion  relations  in  Fig.  1,  only  Che  known  energy  band 
parameters  for  graphite  are  used  with  no  additional  terms  introduced  to  account 
for  the  intercalant-graphite  bounding  layer  interaction.  The  results  for 
Che  p(/3  X  i/3)R30‘*  and  p(2x2)  superlatcices  enable  a  comparison  to  be  made 
between  the  phenomenological  model  and  both  experimental  results  and  theoretica 
first  principles  calculations  for  stage  1  compounds. 

Results  are  shown  in  Fig.  2  for  such  a  comparison  between  the  dispersion 

relations  calculated  for  the  stage  1  empty  Intercalate  layer  model  for 

a  p(»^  X  i/3)R30®  and  p(2x2)R0®  superlattice  and  for  the  first  principles 

calculations  of  Holzwarth  et  al.  [9]  for  C,L1  and  of  Ohno  et  al.  [10] 

o 

for  CgK.  Good  agreement  is  obtained  for  the  bands  related  to  the  graphite 
ir-bands,  showing  that  the  perturbation  of  the  graphite  ir-bands  by  the 
Intercalant  is  small.  To  obtain  agreement  for  the  hybridized  metal  s-bands, 

It  is  necessary  to  fit  the  phenomenological  bands  to  Che  first  principles 
hybridized  s-bands  to  determine  the  overlap  band  parameters  between  the 
graphite  ir-orbltals  and  the  intercalate  orbitals.  One  attractive  feature 
of  the  phenomenological  model  is  that  once  the  band  parameters  are  determined 
for  any  low  stage  compound  (where  the  graphite- Intercalate  interactions  are 


Important) ,  the  dispersion  relations  for  all  stages  can  then  be  evaluated 
without  Introduction  of  additional  band  parameters.  The  phenomenological 
model  thus  provides  a  useful  tool  for  the  calculation  of  the  electronic 
structure  for  high  stage  compounds  where  first  principles  calculations  are 
difficult  to  make  and  the  phenomenological  model  Is  expected  to  be  most 
convergent . 

A  comparison  between  the  dispersion  relations  for  zone  folded  graphite 

(Eq.  4)  and  for  the  empty  Intercalate  layer  model  (Eq.  5)  Is  shown  In  Fig.  3 

for  a  stage  3  compound  for  both  a  p(i/3  x  i/3)R30‘  and  a  p(2x2)  In-plane 

superlattlce.  Independent  of  the  in-plane  super lattice,  the  replacement 

of  a  graphite  layer  by  an  empty  Intercalate  layer  gives  rise  to  only  minor 

differences  In  the  structure  of  the  r-bands,  but  results  In  a  great  reduction 

In  the  k  dispersion  near  the  zero  of  energy  In  Fig.  3.  This  reduction  in 

the  k  dispersion  has  major  consequences  on  the  electronic  and  optical 
z 

properties  of  Intercalation  compounds.  Inclusion  of  a  donor  or  an  acceptor 
layer  does  little  to  perturb  the  general  features  of  the  "empty  lattice"  energy 
band  model  except  for  the  placement  of  the  energy  levels  near  the  zero  of 
energy. 

Finally,  we  show  In  Fig.  4,  dispersion  relations  based  on  the  empty 
Intercalate  layer  model  for  stage  1  to  5  compounds  with  a  p(/3  x  /3)R30® 

In-plane  folding.  On  the  left  we  show  a  plot  of  the  energy  bands  on  an 
expanded  energy  scale  In  the  region  of  the  F  and  A  points  in  the  hexagonal 
Brillouln  zone.  For  this  in-plane  superlattlce,  the  Slonczewski-Weiss-McClure 
bands  fold  into  the  region  of  the  Brillouln  zone  about  FA.  A  comparison 
with  other  In-plane  orderings  shows  that  the  details  of  the  band  structure 
In  the  neighborhood  of  the  Fermi  level  Is  Independent  of  in-plane  order, 
since  no  explicit  Intercalate-graphlte  bounding  layer  Interaction  has  been 
Introduced.  In  this  figure  we  further  note  that  the  dispersion  relations  diffe 
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as  we  move  away  from  the  F  point  In  the  FM  and  FK  directions,  reflecting 
the  trigonal  warping  of  the  energy  bands. 

The  basic  model  for  the  dispersion  relations  presented  In  this  paper  has 
been  applied  to  fit  experimental  Shubnikov  de  Haas  (SdH)  [7]  and  magneto- 
reflection  [8]  data.  In  addition,  structure  observed  in  the  Infrared 
reflectivity  [11,12]  can  be  Identified  with  transitions  between ' energy 
bands  that  are  folded  back  because  of  the  k^  periodicity  due  to  the  staging 
phenomenon.  To  distinguish  between  donor  and  acceptor  compounds  and  to 
obtain  a  quantitative  fit  to  experimental  data.  It  Is  necessary  to  Include 
an  explicit  interaction  between  the  Intercalate  and  graphite  bounding  layers. 

With  regard  to  Shubnikov  de  Haas  measurements,  the  Fermi  surface  cross 
sectional  areas  depend  critically  on  the  position  of  the  Fermi  level.  If 
the  Fermi  level  lies  above  the  highest  rr-band  along  the  HK  axis,  the  SdH 
frequencies  will  group  Into  a  number  of  clusters  equal  to  the  stage  of  the 
compound.  The  splitting  within  the  clusters  is  directly  related  to  the  k^- 
dlsperslon  of  the  energy  levels.  In  the  fit  to  the  SdH  frequencies  observed 
for  stage  5  graphlte-K  [7],  the  Fermi  level  was  placed  so  that  only  4  out  of 
5  possible  electron  pockets  were  occupied,  indicating  a  non-uniform 
z-dependence  of  the  charge  distribution. 

The  z-dependence  of  the  charge  distribution  depends  critically  on  both 
the  position  of  the  Fermi  level  and  the  shift  in  the  bounding  layer  potential 
associated  with  the  intercalate-graphlte  bounding  layer  interaction.  If  the 
ir-bands  are  fully  occupied  (or  fully  empty)  or  exactly  half  occupied  (as  in 
pristine  graphite),  then  the  u-electron  charge  distribution  is  uniform  on 
the  graphite  layers.  Since  each  of  the  energy  bands  corresponds  to  a  different 
contribution  from  each  of  the  real  space  layers,  a  partial  occupation  of  anti- 
bonding  if-bands,  associated  with  raising  the  Fermi  level  through  donor 


8 


intercalation,  gives  rise  to  a  non-uniform  charge  distribution  even  within 
the  empty  intercalate  layer  model.  The  introduction  of  an  intercalate-grnphi tc 
bounding  layer  interaction  further  contributes  to  the  z-dependence  of  the 
charge  distribution. 

The  formalism  developed  for  the  electronic  levels  in  the  Intercalation 
compounds  is  easily  applied  to  the  calculation  of  quantized  Landau  levels  [13]. 
Such  an  analysis  is  needed  for  the  detailed  interpretation  of  the 
magnetoref lection  experiments  [8]  which  show  that  intercalation  gives  rise 
to  small  changes  in  the  graphite  band  parameters.  These  changes  can  be  measured 
quantitatively  over  a  wide  range  of  intercalate  concentrations  and  the  resulting 
band  parameters  can  be  incorporated  Into  the  model  of  Eq,  (5).  The  large 
decrease  found  for  the  k^-dispersion  of  the  electronic  states  is  also  observed 
for  the  Landau  levels,  and  hence  the  magnetic  fleld«>induced  singularities 
in  the  density  of  states  are  enhanced  by  intercalation.  In  addition,  k2-axls 
zone  folding  Introduces  new  singularities  in  the  density  of  states  which  are 
expected  to  give  rise  to  new  series  of  Landau  level  transitions. 

A  great  amount  of  work  yet  remains  to  be  done  to  obtain  electronic  dispersion 
relations  which  can  quantitatively  account  for  the  large  variety  of  detailed 
experimental  information  now  becoming  available  for  the  graphite  inter¬ 
calation  compounds. 
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Table  I  Points  Equivalent  to  the  P  Point  for  Various  In-Plane  Zone  Foldings 


Superlattice 

Equivalent  Points 

Carbon  Atoms 

Per  Unit  Cell 

Notation 

Label 

B  -S-  k 

ro  2  7r  r 

^2 

P(lxl) 

r 

(0,0,0) 

^6 

pCvT  X  /3)R30“ 

r 

(0,0,0) 

S 

P(2x2) 

r 

(0.0,0) 

Ml*  ^2 

^2’  2/3  ’ 

p(/7  X  /7)R19.1'’ 

r 

(0,0,0) 

Wq, 

(7,  0) 

'  7/3 

^18 

P(3x3) 

r 

(0,0,0) 

•^0'  ^1 

0 

0 

p(/l2  X  /l2)R30* 

r 

(0,0,0) 

*0*  •'i 

(3.^.  0) 

(2' 

Tq* 

(3,  0,  0) 

The  equivalent  reciprocal  space 
by  n60*  around  the  k  -axis. 

9. 

vectors  are  given  by  a  rotation  of  Bq 

» u 


FIGURE  CAPTIONS 


Electronic  dispersion  curves  for  ir-bands  in  graphite  and  in 
an  empty  intercalate  layer  model  for  a  stage  1  Intercalation 
compound  with  various  In-plane  zone  foldings,  each  corresponding 
to  the  indicated  commensurate  ordered  Intercalate  superlattice. 

Comparison  between  the  stage  1  empty  intercalate  layer  model  for 

a  p(/3  X  /3)R30“  superlattice  (solid  curves)  and  the  band  structure 

calculation  (dashed  curves)  for  C.Li  by  Holzwarth  et  al.  (Ref.  9). 

o 

The  p(2x2)  structure  (solid  curves)  is  compared  to  the  band 
calculation  (dashed  cur'/es)  for  C.K  by  Ohno  et  al.  (Ref.  10). 

O 

The  first  principles  calculations  (dashed  curves)  also  show 
alkali  metal  s  and  p  bands  and  consider  explicitly  the  aPyo 
Intercalate  stacking  order  in  the  case  of  CqK;  these  features 

O 

are  not  considered  in  the  present  calculation. 

Fig.  3  Electronic  dispersion  relations  for  a  stage  3  compound  assuming 

both  a  p(v?  x  i/3)R30"  and  a  p(2x2)  In-plane  ordering.  For  the 

dispersion  relations  labeled  'folded  graphite’,  k^-axis  zone  folding 

appropriate  to  a  stage  3  compound  has  been  carried  out  while  the 

curves  labeled  'stage  3'  correspond  to  the  replacement  of  a  graphite 

layer  by  an  empty  Intercalate  layer,  thereby  reducing  the  number 

of  energy  levels  and  reducing  the  k  dispersion  of  the  ensrgy  bands. 

z 

Pig.  A  Bnpty  intercalate  layer  7r-band  model  applied  to  a  p(/3  x  /3)R30® 
superlattice  for  stage  1  (  n  ^  5  compounds.  On  the  left  is  shown 
an  expanded  energy  scale  plot  of  the  bands  which  occur  near  the 
Fermi  level  for  each  Intercalation  compound. 


Fig.  1 


Fig.  2 
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Stage  1 


‘  Recent  advances  In  the  preparation  and  characterization 
of  graphite  intercalation  compounds  have  given  rise  to 
detailed  and  accurate  measurements  relevant  to  the  Feral 
surface  topology  [1,2]  and  to  the  electronic  structure  [2,3]. 
Various  authors  [4,5]  have  stressed  that  even  for  low-stage 
compounds  most  of  the  electronic  levels  in  an  energy  range 
of  'X'ilOeV  of  the  Fermi  level  are  directly  related  to  those 
In  pristine  graphite.  Because  of  this  close  relation  to 
pristine  graphite,  the  theoretical  models  which  have  been 
discussed  are  often  expressed  in  terms  of  single  layer 
graphite  bands  [6,7].  The  two-dimensional  models  for  inter¬ 
calated  graphite  [8]  have  been  written  to  include  charge 
transfer  to  the  graphite  and  screening  charge  in  the 
graphite  bounding  layers. 

The  first  principles  calculations  for  the  stage  1  donor 
compounds  C6bi  [9]  and  CgK  [10,11]  have  confirmed  many  of 
the  above  assumptions.  These  calculations  also  emphasize 
that  a  fully  self-consistent  one-electron  band  calculation 
does  not  require  additional  assumptions  about  charge 
transfer  and  screening  charge  density  because  such  quanti¬ 
ties  can  be  calculated  directly  from  the  wave  functions 
used  in  the  calculation.  Nevertheless,  because  of 
practical  considerations  related  to  the  large  number  of 
atoms  in  the  unit  cell  for  the  dilute  compound,  first 
principles  band  calculations  cannot  be  carried  out  for  the 
entire  range  of  stages  and  intercalants  for  which  experimen¬ 
tal  work  has  been  reported.  Therefore  phenomenological 
models,  such  as  the  kj-axls  zone  folded  model  developed  in 
this  paper,  are  needed  to  explain  experimental  data  or  to 
serve  as  an  interpolation  scheme  between  calculated  k 
points  in  first  principles  calculations  and  as  an  extra¬ 
polation  scheme  to  high  stage  compounds. 

The  results  presented  in  this  paper  represent  a 
specialization  of  a  more  general  c-  and  a-axis  zone  folded 
model  presented  elsewhere  [12,13].  The  general  model  shows 
that  the  form  of  the  levels  near  the  Fermi  surface  reduces 
to  the  Slonczewski-Weiss-NcClure  (SWIcC)  form,  independent 
of  the  In-plane  order.  The  SWIcC  effective  Hamiltonian 
valid  for  kg  on  or  near  the  HKH  axis  is  written  [14,13]  as 


«1 

0 

«13 

* 

«13 

* 

0 

* 

=2 

* 

«23 

-«23 

o 

«13 

«23 

®3 

* 

”33 

"l3 

””23 

«33 

=3 

In  which  the  tc_  dependence  of 

the 

matrix 

the  SWilcC  Hamiltonian 

[14]  pr 

the 

full  3- 

(1) 


expansion  [16]  of  the  interacting  x-bands^ 

The  SWMcC  Hamiltonian  is  written  [14]  in  a  represen¬ 
tation  whose  symmetrized  Bloch  functions  are  given  by 
Yjj  -  (a+a')//2,  -  (a-a’)//2.  -  b.  -  b'. 
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These  functions  transform  as  the  irreducible_^representatlons 
of  the  group  of  the  Ics  wave  vector  at  point  ks  and  are 

a(ka)  -  —  5  e  ^*"8  •  d  rp  (r  -  d) 

^  /S  d  * 

a'(tt  )-  —  J  e  il>  (r  -  d  -  t  ,) 

*  d  *  ^  (2) 


b(k )  -  ^  5 

*  i/if  d 


H  z  o 


(k  )  -  —  ?  e  -  d  -  t.,) 


>diere  the  origin  is  at  atom  A  in  the  A  plane  and  ty^*,  tg* 
tg*  are  vectors  from  the  origin  to  the  other  atoms  in  the 
graphite  unit  cell.  The  Intercalation  compound  is  modelled 
using  a  kj-axis  zone  folding  technique.  The  Hamiltonian  is 
%nrltten  in  terms  of  I  pairs  of  graphite  planes  ABAB...AB, 
where  the  basic  repeat  distance  Is  determined  by  the  number 
of  planes  in  the  c-axls  repeat  distance  of  the  Intercalation 
compound,  considering  for  the  moment  the  intercalate  layer 
to  be  Indistinguishable  from  graphite.  The  zone-folded 
Hamiltonian  is  taken  to  be 


Hr  )• 

folded  s 


H  (k  )  0, 

O  8 

o  . 


_  .  f-1  TT 

H  (k  +  '-r—  ^  z) 
0  8  I  c  ' 

O 


where  each  of  the  blocks  is  the  4xA  SWicC  Hamiltonian  given 
by  Eq.  (1)  and  J!,=n+1  for  even  stage  or  (n+l)/2  for  odd  stage 
compounds,  where  n  is  the  stage  index.  The  basis  functions 
for  the  folded  Hamiltonian,  Eq.  (3),  are  given  by  Eq  •  (2)^ 
with  the  appropriate  translations  in  T5-space,  and  Co®3.35A. 

In  order  to  consider  explicitly  the  effect  of  inter¬ 
calation,  the  Hamiltonian  given  by  Eq.  (3)  must  be  trans¬ 
formed  to  a  layer  representation  which  is  written 
explicitly  by 

■  "  "tolded'*,’  ») 

where  the  unitary  transformation  is  derived  from  the  basis 
functions  given  in  Eq.  (2),  yielding 


In  which  the  blocks  ,...are  2x2  matrices 

having  layer  subscripts.  ^  ^ 

The  unitary  transformation  of  does  not  change 

the  eigenvalues  of  the  Hamiltonian.  Thus  the  energy  bands 
from  the  folded  SWMcC  Hamiltonian  in  the  layer  representation 
of  Eq.  (6)  are  identical  to  the  SWMcC  energy  bands.  However, 
in  the  layer  representation,  the  effect  of  intercalation  is 
easily  included  by  substitution  of  matrix  blocks  for  the 
intercalate  species  in  place  of  the  appropriate  graphite 


layer.  For  the 

Hamiltonian  is 
• 

odd  stage 

n  compounds. 

the  resulting 

“xx 
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“gic 
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where  X  is  written  for  (n-l)/2,  and  X  denotes  an  intercalate 
layer.  For  the  even  stage  compounds  the  effective 
Hamiltonian  is 


V 

“v  H 


I  X  I 

where  p-n/2  and  v-(n/2)+l.  If  the  stacking  of  the  inter¬ 
calate  layers  is  different  in  successive  layers  (e.g.  a,Br 
Che  Hamiltonian  can  be  generalized  by  considering 
explicitly  a  larger  unit  cell  and  appropriate  interaction 
teru  Co  apliC  the  zone-folded  levels. 


In  this  paper  we  consider  the  impllca'tions  of  the  forn 
of  Eq.  (8a)  and  (8b)  by  an  explicit  dlagonalizatlon  of  these 
matrices,  first  for  an  "empty"  lattice  in  which  the  inter¬ 
calate  layer  is  "vacuum"  (l.e.  Hxx  **  HxAi  "  HxBj  =  0)  and 
then  for  a  generalized  donor  or  acceptor  compound  in  which  a 
dispersionless  intercalate  layer  interacts  with  the  graphite 
bounding  layer.  Because  of  the  non-equivalent  bounding  layer 
carbon  sites,  the  Fourier  expansion  Haniltonian  for  the 
Intercalation  compound  includes  a  constant  energy  shift  for 
the  bounding  graphite  layer,  Po*  well  as  k-dependent 
Interaction  terms.  The  results  for  a  single  parameter  model 
are  shown  in  Fig.  1.  The  case  of  the  "empty  lattice" 

(y^  ■  0)  Is  given  in  the  center  of  each  set.  Results  for  a 
bounding  layer  energy  shift  of  Po  0.3eV  are  presented  on 
the  left  side  labeled  acceptors,  and  for  Pq  =  -0.3eV  on  the 
right  side  labeled  donors.  The  position  of  the  intercalate 
level  Ex  Is  assumed  to  lie  below  the  region  of  the  SWMcC 
ir-bands  for  acceptors  upshlfting  the  bounding  layer  energy 
levels.  Correspondingly  is  above  the  ‘n'-bands  for  donors 
giving  a  downshift  to  the  bounding  layer  levels.  The 
interaction  of  0.3eV  was  suggested  by  the  preliminary 
analysis  of  Shubnikov-de  Haas  data  [1,2].  The  plots  in  Fig. 

1  are  made  using  the  SWMcC  band  parameters  for  pristine 
graphite  [17]  given  in  Table  1  and  a  single  interaction 
parameter  p©*  The  magnitude  of  Pq  is  expected  to  differ  for 
each  Intercalant,  but  provided  the  in-plane  intercalate 
density  Is  independent  of  stage,  Pq  will  also  be  independent 
of  stage.  Thus,  the  dispersion  curves  given  in  Fig.  1  for 
stages  1  ^  n  ^  8  depend  on  the  introduction  of  a  single 
parameter  Pq. 

The  withdrawal  of  a  graphite  layer  greatly  reduces  the 
dispersion  of  the  bands  in  the  kz  direction.  Additional  kz 
dispersion  can  be  introduced  by  inclusion  of  kz-dependent 
terms  in  the  Intercalate-graphite  bounding  layer  interaction. 
For  stages  n  >  2,  the  near  degeneracy  of  the  E3  graphite 
levels  shown  for  the  "empty  lattice"  bands  in  Fig.  1  is 
partially  lifted  by  Pq  through  non-rigid  band  shifts,  though 
a  number  of  the  E3  levels  still  cluster  and  retain  their 
near  degeneracy.  However,  for  stages  n»l  and  2,  the  effect 
of  Po  is  merely  to  produce  a  constant  energy  level  shift, 
which  is  physically  reasonable  since  in  these  cases  there  are 
only  graphite  bounding  layers.  Note  that  for  higher  stage 
compounds  (n  ^  4)  in  Fig.  1,  both  electron  and  hole  pockets 
can  be  formed  if  the  interaction  energy  is  comparable  with 
the  Fermi  energy  shift. 

Landau  level  transitions  between  the  nearly  degenerate 
E3  levels  a''e  observed  in  the  magnetoref lection  spectra  [3] 
which  are  well  fit  by  the  SWMcC  model  with  small  but  measur¬ 
able  changes  in  the  SWMcC  band  parameters  upon  intercalation. 
These  modified  values  of  the  band  parameters  can  then  be 
used  in  applying  F.q.s.  (da)  and  (8b)  to  the  interpretation  of 
experimental  data  such  as  the  de  Haas- van  Alphen  effect. 
Though  the  parameter  changes  result  in  small  perturbations 
of  Che  SWMcC  bands,  near  the  HK  axis,  larger  interactions 


could  occur  in  other' regions  of  the  Brillouin  zone,  e.g.  the 
M  point.  The  energy  bands  of  Fig.  1  show  extremal  cross 
sectional  areas  at  the  H  and  K  points  [2]  which  can  be 
calculated  and  compared  with  measured  DHVA  frequencies  and 
effective  masses.  The  band  model  of  Eq.  (8a, 8b)  is  also 
amenable  to  Landau  level  calculations  using  either  raising 
and  lowering  operator  techniques  [18]  or  semi-classical  Bohr- 
Sommerfeld  quantization  [19],  both  of  which  have  been  applied 
to  graphite. 

The  bands  shown  in  Fig.  1  are  all  graphite  ir-bands. 

The  levels  represented  by  the  intercalate  Hamiltonian  Hxx 
are  not  Included  in  this  presentation  of  results  in  order 
to  show  features  of  the  graphite  tr-bands  which  are  common 
to  all  Intercalants.  The  intercalate  block  HxX  Is  eliminated 
from  Eqs.  (8a)  and  (8b)  using  the  LBwdin  summation.  This 
gives  the  addition  of  a  term  -(Hxa^Hxa- ) /(Hyx  ~  £)  to  the 
larger  block  in  Eq.  (8)  where  e  is  the^eigenvalue.  The  sign 
of  this  contribution  is  negative  for  donors  when  the  inter¬ 
calate  layers  are  above  e  and  positive  for  acceptors  where 
the  Intercalate  layers  are  below  e.  The  Fermi  level  can  be 
calculated  by  assuming  explicitly  the  fractional  charge 
transfer  per  Intercalate  atom  or  molecule.  Incomplete 
charge  transfer  to  the  SWMcC  bands  implies  that  the  inter¬ 
calate  bands  are  occupied  or  that  in-plane  zone-folding  has 
resulted  in  the  occupation  of  Tr-bands  at  other  points  in 
the  zone  such  as  the  M-point  [13].  These  additional  levels 
should  then  be  Included  explicitly  in  Eqs,  (8a, 8b). 

The  electronic  states  associated  with  the  intercalate 
layer  can  be  explicitly  incorporated  in  the  Hxx  block  of  the 
Hamiltonian,  given  in  Eqs.  (8a, 8b).  The  alkali  metal 

donor  compounds  require  Levels  associated  with  s  and  p 
orbitals,  whereas  the  acceptor  compounds  generally  require 
Che  use  of  more  localized  electronic  states.  These  consider¬ 
ations  for  the  electronic  levels  in  the  pristine  intercalate 
material  suggest  that  donor  compounds  should  exhibit  more  kz 
dependence  in  the  coupling  blocks,  HxBq,  and  more  k^  disper¬ 
sion  Chan  Che  acceptor  compounds.  The  observed  conductivity 
anisotropies  in  donors  and  acceptors  is  consistent  with  this 
basic  different  in  the  intercalate  electronic  levels. 

This  simple  extension  of  the  phenomenological  graphite 
Hamiltonian  features  the  close  relation  bewteen  pristine 
graphite  and  the  well-staged  intercalation  compounds.  The 
model  is  based  on  a  kz-axis  Brillouin  zone-folding  and  makes 
no  use  of  a-plane  superlattice  periodicity.  The  same 
Hamiltonian  with  no  additional  parameters  can  be  used  Co 
model  all  stages  (plots  are  shown  in  Fig,  1  for  stages 
1  <  n  <  8).  The  model  can  be  used  for  the  standard  intro¬ 
duction  of  Landau  level  quantization  and  thus  provides  a 
direct  way  of  applying  quantum  oscillatory  phenomena  to  a 
quantitative  determination  of  the  electronic  structure  of 
Intercalation  compounds. 
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Figure  Captions 


Fig.  1  E(k}  for  graphite  intercalation  compounds  along  the 
the  edge  of  the  hexagonal  Brlllouin  zone  where  the 
K-polnt  corresponds  to  k2*0.  The  E(^)  for  the  odd- 
stage  compounds  are  plotted  in  (a)  for  stages 
n“l,3,5,7  where  the  H-point  corresponds  to  tt/Ic 
where  Ic  is  the  repeat  distance  in  the  c-directlon. 
The  E(k)  for  the  even  stages  is  plotted  in  (b)  for 
n“2,4,6,8  where  the  H-point  corresponds  to  n/ZIc- 
An  Intercalate  graphite  bounding  layer  interaction 
appropriate  to  donors  is  shown  on  the  right  hand 
side  of  the  figure  and  that  appropriate  to 
acceptors  is  shown  on  the  left.  For  each  stage, 
the- center  panel  is  for  the  "empty  lattice"  with 
no  graphite- inter calant  interaction.  Where 
degeneracies  (or  near  degeneracies)  occur,  the 
number  of  levels  is  indicated.  For  the  even  stages 
the  levels  pair  away  from  the  HK  axis  except  for 
the  double  pairings  labeled  4. 
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The  kz^axis  zone  folded  Hamiltonian  based  on  pristine 
graphite  energy  bands  with  an  empty  intercalate  layer 
accounts  for  both  observations  of  the  Fermi  surf^ice 
made  by  the  Shubnikov-de  Haas  technique,  wliich  s!iow 
large  changes  with  intercalation,  and  of  the  magnetic 
energy  level  structure  by  the  magnetoref lection 
technique,  where  very  small  changes  are  found.  This 
model  has  implications  for  transport  and  optical 
properties  for  donor  intercalation  com.pounds. 

Because  of  the  unusual  properties  of  graphite  intercalation 
compounds  and  their  close  connection  to  the  semiconductor  super¬ 
lattice  structures  now  being  fabricated  by  molecular  beam  epitaxy, 
graphite  intercalation  compounds  have  attracted  attention  in  the 
semiconductor  community.  Analogous  to  the  transport  behavior 
observed  for  the  semiconductor  surjorlattice  structures,  the  iiigh 
electrical  conductivity  of  graphite  intercalation  compounds  is 
achieved  by  charge  transfer  from  the  intercalate  layers  where  the 
mobility  is  low  to  the  graphite  layers  \>/hore  the  mobility  is  very 
high.  The  superlattice  structure  created  by  molecular  beam  epitaxy 
growth  is  achieved  in  graphite  intercalation  compounds  by  the  staging 
mechanism  which  results  in  intercalated  graphite  with  sequential 
intercalant  layers  separated  by  n  graphite  layers  where  n  denotes 
the  stage.  For  intercalated  graphite  the  staging  periodicity  is  long 
range  and  is  a  dominant  symmetry. 

The  electronic  structure  is  determined  by  the  basic  symmetry  of 
graphite  as  modified  by  the  staging  periodicity,.  This  structure  is 
here  studied  by  the  Shubnikov-de  Haas  effect  for  a  determination  of 
the  Fermi  surface  and  by  the  magnetcref lection  technique  which  ..s 
sensitive  to  the  quantized  electronic  energy  levels  in  a  magnetic 
field.  A  phenomenological  model  based  on  crystal  syit'.metry  is  used  to 
interpret  and  to  interrelate  the  experimental  observations  made  by 
the  two  techniques.  Application  of  the  calculational  teclmiques 
level cped  for  the  inter  lated  graphite  superlattice  structures  to  the 
conductor  sunc’' latt  ice  structures  coul  I  be  of  significance. 

;on  of  only  small  changes  in  the  electronic  band  para- 
»  r*-  ,r-  •  I  nt  er :  :i  1  at  ton  suggests  a  theory  for  the  electronic 
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Fig.  1  Electronic  band  structure 
in  the  vicinity  of  the  HK  axis  for 
a  stage  5  "empty  intercalate  layer" 
model  for  a  graphite  intercalation 
compound.  The  Fermi  energy  is 
determined  to  fit  the  observed  SdH 
frequencies.  For  the  indicated 
Fermi  leveli  four  of  the  conduction 
bands  are  partially  occupied  and 
give  rise  to  Fermi  surfaces  and 
carrier  pockets. 


Fig.  2  Magnetoreflection  spectra 
using  (+)  circular  polarization 
for  an  acceptor  (FeCl3  sta^e  ?} 
and  for  a  donor  compound  (P.b 
stage  6)  at  the  indicated  photon 
energies.  For  comparison,  traces 
for  graphite  are  shown.  The 
resonances  are  specif.i.eu  b)’ 
quantum  numbers  for  the  initial 
and  final  states  (Ref.  3), 


interaction  is  neglected  [2],  These  results,  shown  for  the  3 
conduction  and  5  valence  bands  near  the  H>K  zone  edge,  demonstrate 
that  intercalation  greatly  decreases  the  kz  dispersion.  The  place¬ 
ment  of  the  Fermi  level  in  Fig.  1  results  in  a  partial  occupation 
of  the  TT-bands  which  gives  rise  to  a  z-dependen't  charge  distribution. 

Detailed  information  on  the  form  of  these  energy  bands  is 
provided  by  the  mrgnetoreflcction  experimjnt  which  probe?  resonant 
transitions  from  occupied  magnetic  energy  levels  below  the  Fermi 
level  Ep  to  unoccupied  levels  above  Ep.  Measurements  have  been  made 
on  well-staged  donor  (stages  n^6)  and  acceptor  (stage  n?4)  compounds 
[3],  characterized  by  (OOt)  dif fractograms  to  establish  their  stage- 
index  and  staging  fidelity.  Resonant  structures  observed  in  the 
magnetoreflection  spectra  are  identified  with  Landau  level  transi¬ 
tions  between  magnetic  energy  levels  associated  with  the  H3  bands 
which  in  graphite  are  degenerate  at  the  K-point,  The  close  similnrit 
of  the  magnetoreflection  spectra  in  pristine  graphite  and  in  those 
donor  and  acceptor  comounds  indicates  tliat  the  eJectronic  strticnue 
of  the  n -bands  near  Ep  remains  highly  graphitic  upon  i ntcrca ia t i on . 


Table  1  Results  of  analysis  of 
magnigtoref lection  experiments 


Inter¬ 

calant 

Stage 

A 

m 

c 

A 

m 

V 

(HOPG) 

00 

0.056 

0.084 

25.1 

AlCl, 

8 

0.056 

0.076 

26.0 

6 

0.054 

0.076 

26.6 

FeCl. 

7 

0.055 

0.079 

26.1 

5 

0.054 

0.075 

26.6 

Rb 

6 

0.045 

0.065 

31.7 

separation.  A  summary  of  typical 
strong  support  for  the  close  relat 
the  intercalation  compounds  and  in 


Analysis  of  the  observed  spectra 
on  the  basis  of  a  two-band  mod'.-l 
for  the  K-pcint  valence  and 
conduction  bands  yields  tlic 
valence  and  conduction  band 
effective  masses  iiiv  and  Since 

the  observed  magnetoreflection 
spectra  in  the  intercalation  com¬ 
pounds  are  well  satisfied  by  the 
form  of  the  Slonczewski-Weiss - 
McClure  (SV.'McC)  energy  level  model 
for  graphite,  analysis  of  these 
data  yields  the  St\fMcC  band  para¬ 
meter  combination  (yq/Yi )  which 
determines  the  Landau  levei 
results  is  given  in  Table  1,  giving 
ion  between  the  graphite  it -bands  in 
pristine  graphite. 


To  study  the  Fermi  surface,  Shubnikov-de  Haas  (SdH)  measurements 
of  the  transverse  magnetoresistance  have  been  made.  Particular 
emphasis  was  given  to  measurements  on  encapsulated,  well-staged  alkali 
metal  donor  compounds  [2] .  The  four-point  method  was  used  to  study 
the  a-plane  transverse  magnetoresistance  in  the  temperature  range 
1.4<T<4.2K  and  in  magnetic  fields  up  to  15  Tesla.  Data  acquisition 
was  by  computer,  and  the  data  were  manipulated  to  obtain  a  Fourier 
power  spectrum  of  resistance  vs  1/H,  thereby  yielding  the  frequoicies 
of  SdH  oscillations,  which  are  related  to  the  extremal  cross  sections 
of  the  Fermi  surface  [2],  The  Fourier  power  spectra  for  the  SdH 
frequencies  for  potassium  stages  2,5  and  8  are  shown  in  Fig.  3.  Of 
particular  significance  is  the  stage  dependence  of  the  observed 
spectra,  in  agreement  with  results  by  Suematsu  et  al.  [4] .  A  detailed 
study  was  made  for  several  stage  5  samples  and  the. observed  SdH 
frequencies  compared  with  predictions  of  the  empty  intercalate  layer 
model  of  Fig,  1,  with  the  Fermi  level  placed  to  fit  the  Sdll  frequen¬ 
cies.  It  should  be  noted  that  the  empty  intercalate  layer  model  uses 
only  graphite  band  parameters,  which  have  been  determined  previously. 
Thus  the  good  agreement  obtained  between  the  calculated  and. observed 
SdH  frequencies  in  Table  2  gives  strong  support  that  the  energy  levels 
close  to  Ef  are  very  similar  to  the  basic  ir-bands  of  pristine  grapliitc 
in  agreement  with  the  magnetoreflection  results  discussed  above.  The 
large  change  in  the  SdH  frequencies  and  their  stage  dependence  is  due 
to  the  stage  dependence  of  the  number  of  conduction  and  valence  bands, 
of  the  k2*axis  zone-folding  procedure  and  of  the  position  of  the  ret:'  i 
level  relative  to  the  K-point  band  edge.  It  is  of  interest  to  note 
that  light  masses  similar  to  those  in  pristine  graphite  are  found  for 
the  small  carrier  pockets,  with  significantly  larger  masses  for  the 
larger  carrier  pockets  which  probably  contribute  more  significantly 
to  transport  properties.  It  is  furthermore  likely  that  other  carrier 
pockets  of  lower  mobility  occur  elsewhere  in  the  zone,  tliereby 
accounting  for  the  proper  charge  transfer. 

From  the  Shubnikov-de  Haas  results  and  the  model  for  the  energy 
bands  we  conclude  that  the  transport  properties  of  graphite  inter¬ 
calation  compounds  must  be  interpreted  in  terms  of  a  multi -pocket 
model,  with  contributions  from  each  of  the  occupied  carrier  pockets 
(see  Fig.  1  and  Table  2).  Explicit  use  of  an  energy  band  model  which 
can  be  applied  to  any  stage  and  intercalant  provides  a  basis  for 
calculations  of  the  transjiort  properties  for  intercalated  gvnp’iLte. 


Fig.  3  Shubnikov'de 
Fourier  transform  pouer 
spectra  for  stages  2,5  and  6 
graphite-K.  These  power 
spectra  were  obtained  by  a 
Fourier  transform  of  an 
experimental  resistivity  vs 
1/H  trace  for  magnetic  fields 
H<15  Tesla.  The  peaks  in  the 
power  spectra  correspond  to 
SdH  frequencies,  whicli  are 
given  in  Tesla  and  the  saine 
scale  is  used  for  eacii  stage. 


Table  2  Fermi 

surface 

parameters 

associated  with 

stage  5  graphite-K 

K-point  Band 

Designation 

Fermi  Surface 

^1 

Parameters 

m  /m^ 

0.143 

0.115 

0.0828 

0.0511 

n^(xl0  cm  ) 

2.36 

1.76 

0.975 

0.138 

SdH  Frequencies 
Calculated 

(Tesla) 

401 

300 

163 

26.7 

> 

Observed  j 

453 

290 

191 

24 

(Tesla) 

430 

267 

152 

13 

243 

135 
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Electronic  Energy  Bands  and  Fermi  Surfaces  in 
Graphite  Intercalation  Compounds. t  G.  DRESSELHAUS,* 

S.Y.  LEUNG' and  M.  SHAYEGAN,  M.I.T.— A  model  for  the 
electronic  energy  bands  for  graphite  intercalation 
compounds  is  presented,  based  on  the  3-dlmenslonal 
Fourier  expansion  for  the  fr-bands  in  pristine  graphite. 
The  model  focuses  both  on  staging  periodicity  which  is 
created  by  zone  folding  in  the  direction  and  an  in¬ 
plane  superlattice  which  introduces  an  in-plane  zone¬ 
folding.  The  effect  of  Intercalation  is  treated  in 
the  Hamiltonian  by  the  removal  of  a  graphite  layer  and 
Its  replacement  by  an  intercalate  layer,  resulting  in 
a  greatly  reduced  dispersion  with  only  minor  changes 
in  the  In-plane  dispersion.  Graphite-lntercalate 
interactions  are  introduced  to  yield  the  appropriate 
kz  dispersion  indicated  by  experiment.  This  model  is 
used  to  calculate  extremal  areas  of  the  Fermi  surface 
which  are  compared  with  SDH  and  DHVA  experiments  on 
low  stage  alkali  metal  compounds. 
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Angular  Dependence  of  SdH  Frequencies  in 
Alkali-Metal  Graphite  Intercalation  Compounds. t* 

M.  SHAYEGAiN  and  G.  DRESSELHAUS'r,  M.I.T.—  The  different 
stage  dependence  of  the  Fermi  surface  reported  for 
donor  and  acceptor  compounds  has  been  a  matter  of 
Interest  and  controversy.  We  report  the  angular 
dependence  (d=angle  between  magnetic  field  and  c-axis) 
for  Shubnlkov-de  Haas  (SdH)  frequencies  ^0 (Tesla)  in 
veil-staged  (stages  n^4)  graphite-K  and  graphite-Rb 
compounds.  The  observed  Oq  for  potassium  stage  2  are 
135T,  150T,  280T  and  310T,  with  the  main  at  150T. 

The  6  dependence  of  these  is  consistent  with  a 
cylindrical  Fermi  surface  along  the  c-axls.  Rubidium 
samples  show  significantly  different  features. 

Rubidium  stage-2  has  its  main  P.q  at  28T  with  two  fio  of 
lower  intensity  at  220T  and  260T.  The  0  dependence 
for  these  fi©  is  not  compatible  with  cylindrical  Fermi 
surfaces.  Comparison  is  made  to  DHVA  frequencies  for 
graphite-K  compounds  reported  by  Tanuma  et  al.  The 
results  for  K  and  Rb  donor  intercalants  Is  contrasted 
and  compared  with  observations  in  acceptor  compounds. 
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ABSTRACT 

Raman  spectroscopy  has  been  used  to  study  Br-Br  bond  stretch  modes, 
both  for  bromine  adsorbed  on  Grafoil  between  -120  and  -ISS'^C,  and  for 
bromine  in  the  gas  phase  above  a  room  temperature  pyrolytic  graphite 
sample.  The  low  temperature  spectra  include  broad  (FIVHM— 24  cm-1)  peaks 
at  269  and  277  carl,  and  a  narrower  peak  (  —  4  cm-1)  at  274  cm-1.  These 
three  peaks,  which  are  seen  individually  at  different  temperatiires  ^md 
bromine  pressures,  are  tentatively  correlated  with  the  pressure/temperature 
adsorbed  phase  diagram  from  Lander  and  Morrison's  LEED  study.  The  data  do 
not  cover  a  sufficient  ratnge  of  pressures  and  temperatures  to  make  a 
positive  identification  of  the  observed  modes  with  the  adsorbed  phases. 

At  excitation  wavelengths  shorter  than  5100  Angstroms,  the  room 
temperature  Raman  spectra  are  featureless  at  all  bromine  pressures  studied. 
In  particvilar,  the  expected  Br-Br  stretch  modes  corresponding  to  adsorbed 
Br2  molecules  at  gas  pressures  near  the  bromine-graphite  intercalation 
threshold  (18  Torr)  are  not  observed. 

At  excitation  wavelengths  longer  than  5100  Angstroms,  sharply 
resoneuit  room  temperature  Raman  spectra  are  observed.  Relatively  simple 
spectra  excited  by  argon  and  krypton  laser  lines  are  identified  as  arising 
when  the  photon  energy  matches  the  energy  separation  between  a  level  of 
the  ground  state  and  a  level  of  the  ^TTcu'^  excited  state.  Dye  laser 

spectra  over  the  range  5722  A  to  6360  X,  which  are  much  more  complex,  are 
nevertheless  shown  to  arise  from  the  same  source. 

A  background  chapter  surveys  the  literature  concerning  bromine  inter¬ 
calation  and  adsorption  on  graphite,  as  well  as  resonant  scattering  from 
rotational  energy  levels  of  Br2  in  the  gas  phase. 

An  experimental  chapter  discusses  the  equipment  used  in  the  light 
scattering  experiments. 

An  experiment  is  proposed  to  study  fluorescence  from  resonant 
absorptions,  in  order  to  deduce  tlte  energies  of  high  rotational-vibrational 
levels  of  the  ground  state. 
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1.  Introduction 


Graphite  intercalation  compounds  represent  a  class  of  compounds  with  many 
similarities  to  the  superlattice  semiconductors  prepared  by  molecular  beam 
epitaxy.  Intercalated  graphite  has  recently  attracted  considerable  atten¬ 
tion  because  the  addition  of  an  intercalant  can  dramatically  change  the  elec¬ 
tronic  properties  of  the  host  graphite  material.  For  example,  with  the  ad¬ 
dition  of  AsF^,  an  intercalation  compound  with  room  temperature  conductivity 
comparable  to  that  of  copper  can  be  achieved  [1],  while  the  addition  of  an 
alkali  metal  can  result  in  a  superconductor  [2].  Similar  to  the  behavior  in 
modulation-doped  superlattices,  high  conductivity  in  intercalated  graphite 
is  achieved  by  a  charge  transfer  from  the  intercalate  layer,  where  the  mo¬ 
bility  is  low,  to  the  high  mobility  graphite  layers. 

High  magnetic  fields  allow  use  of  the  magnetoreflection  technique  to  ob¬ 
tain  detailed  information  on  the  electronic  dispersion  relations  x^ithin  a 
few  hundred  meV  of  the  Fermi  level  end  use  of  the  Shubnikov-de  Haas 

effect  and  related  quantum  oscillatory  phenomena  to  provide  information  on 
the  Fermi  surface  [5-8].  Initially  it  seemed  difficult  to  reconcile  the  re¬ 
sults  obtained  by  these  two  high  field  techniques  because  the  magnetoreflec¬ 
tion  experiment  provided  strong  evidence  that  the  electronic  structure  near 
the  Fermi  level  remained  highly  graphitic  upon  intercalation,  x^hile  the 
Shubnikov-de  Haas  experiment  indicated  that  large  changes  in  the  Fermi  sur¬ 
face  occur  through  intercalation.  Using  a  recently-developed  phenomenologi¬ 
cal  model  for  the  electronic  structure  of  intercalated  graphite  [9],  consis¬ 
tency  between  the  magnetoref lection  and  Shubnikov-de  Haas  results  is  demon¬ 
strated.  This  model  also  provides  a  basis  for  the  calculation  of  the  mag¬ 
netic  energy  level  structure  which  is  needed  to  explain  the  magnetoref lection 
spectra  in  detail. 

Graphite  intercalation  compounds  are  prepared  by  the  insertion  of  atomic 
or  molecular  layers  of  a  different  chemical  species  between  layers  of  the 
graphite  host  material  [10].  In  this  insertion  or  intercalation  process,  the 
layer  structure  of  graphite  Is  almost  completely  preserved  because  of  the 
acrong  intrapl.nnar  binding  in  graphite,  and  likewise  the  intercalate  layer  is 
closely  related  to  a  layer  in  the  parent  solid.  Structurally,  graphite  in- 
tcrcalatlnn  compounds  exhibit  tlxe  rt-iii.'irkable  phenomenon  called  «  w'lere 

by  the  Intercalate  layers  .irc  periodically  arr.-inged  in  a  matrix  of  gr.tphiie 
layers  and  are  rhnr.xetor  1 /ed  hv  a  Index  n,  d.  nut  Ing  the  numhe-  i-  :  t  i 
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S5unerl.it cicc  structure  similar  to  chat  prepared  in  simiiconducf ors  t'so 
liioiecular  beam  epitaxy  teclmique.  Staging  is  widespread  among  granhice  in¬ 
tercalation  compounds  and  well-staged  siimnles  can  be  prepared  up  to  hig’.i 
stage  (n'vlO).  The  stage  determination  and  the  c-axis  repeat  distance  T  .  .nre 
obtained  from  (OOi)  x-ray  di!  fraccograms.  Because  of  the  strong  st.jge  de  -ea- 
dence  found  for  the  Shubnlkov-de  Ha.as  (SdH)  frequencies  in  the  donor 
calation  compounds  in  this  work,  it  is  necessary  to  carry  out  tlie  '.■..easur o- 
ments  on  single-staged  samples.  By  careful  preparation  of  well-staged  co'’.,- 
pounds  it  is  possible  to  achieve  the  resonance  condition  Wj,t>1  reuu''red  for 
the  observation  of  Landau  level  phenomena.  For  both  the  magnetoreflection 
and  Shubnikov-de  Haas  studies,  the  samples  were  prepared  by  the  tv.'o-zone 
growth  technique  [11]  and  the  stage  index  was  determined  from  (OOi)  x-ray 
dlffractograms. 


2.  Magnetoreflection  Experiments 


The  magnetoreflection  measurements  were  made  at  nearly  normal  incidence  to 
the  c-face  of  the  intercalated  graphite  samples  using  the  Faraday  geometry 
with  magnetic  fields  in  the  range  0<H£l5  Tesla  [4].  Magnetoref lection 
traces  were  taken  at  constant  photon  energy  in  the  range  0.1  eV  <  Trj  <  0.5eV, 
using  circularly  polarized  light  obtained  with  a  gold-xoire  grid  and  a  Csl 
Fresnel  rhomb.  Measurements  were  made  using  a  cold  finger  dewar  operating 
at  4.2K. 


The  observed  magnetoreflection  spectra  are  in  most  cases  qualitatively 
similar  to  those  of  pristine  graphite,  but  show  differences  with  regard  to 
resonant  magnetic  fields  and  in  some  cases  also  differences  in  resonant  line- 
shapes.  Illustrative  spectra  are  sho\.m  in  Fig.  1  for  an  acceptor  compound 


Graphite 
hw  0.338  eV 


Rb  Sioge  6 
t»w  0.335  eV 


Fig.  1  Magnetoreflcctlon  spectrum 
using  (+)  circular  polarization  for 
an  acceptor  compound  (FeCl-j  stage  7) 
at  a  photon  energy  ■hu>  =  0.295  eV  and 
for  a  donor  compound  (Kb  stage  6)  at 
hu*  0.335  eV.  For  comparison,  traces 
for  graphite  arc  shown  at  comparable- 
pltoLon  energies.  The  resonances  ari 
specified  by  the  qu.intum  numbers  for 
the  initial  and  fin.Tl  st.itcs. 
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(FeCl3  sta^je  7)  and  a  donor  conipound  (Rb  stage  6).  For  comparison,  spocLia 
are  also  included  lor  pure  grapuile  at  similar  plioton  energies.  Kacli  struc¬ 
ture  In  Fig.  1  is  identified  with  a  K-point  Landau  level  interband  transition 
specified  by  the  quantum  numbers  for  the  initial  and  final  states.  For  ail 
compounds  that  have  been  studied,  tl\e  resonant  magnetic  field  for  a  riven 
transition  shifts  to  increasingly  lower  fields  with  increasing  intercaiate 
concentration  for  both  (+)  and  (-)  senses  of  circular  polarization.  The 
shifts  of  the  resonant  magnetic  fields  for  a  given  transition  are  much  larg¬ 
er  for  donor  compounds  than  for  acceptors  of  comparable  stage.  The  inten¬ 
sity  of  the  resonances  (magnitude  of  the  percent  reflectivity  change)  is  re¬ 
duced  with  respect  to  tliat  in  graphite  because  of  a  reduction  in  relaxation 
time  and  an  increase  in  tlie  magnitude  of  the  zero  field  reflectivity  !U0)  . 

For  example,  at  a  photon  energy  of  ■htjj'uO.2  eV,  R(0)  is  found  to  be  't0,7, 

0.8,  0.9  respectively  for  pristine  graphite,  graphite-FeCl3  stage  7,  and 
graphite-Rb  stage  6. 

Resonant  Landau  level  transitions  are  observed  and  are  analyzed  by  taking 
spectra  over  a  wide  range  of  photon  energies.  A  suiranary  of  such  results 
for  a  typical  intercalation  compound  is  shown  in  Fig.  2,  in  which  a  compari¬ 
son  with  results  for  pristine  graphite  is  included.  From  this  analysis  it  is 


Fig.  2  Summary  of  resonant  magnetic 
fields  ’’fan  chart"  for  photon  ener¬ 
gies  in  the  range  0. 110-^uj<0 . 520  eV 
for  a  FeCl3  stage  7  sample  upr:;g  (+) 
circularly  polarized  radio t.' on.  ’’'o 
emphasize  the  similarities  but  mea¬ 
surable  differences,  a  comparison  of 
the  results  for  this  compound  and 
for  graphite  is  presented  in  the 
same  photon  energy  range. 

shown  that  for  graphite  intercalation  compounds  with  stages  the  elec¬ 

tronic  structure  within  a  few  hundred  millivolts  of  the  Fermi  level  is  well 
described  by  the  graphite  Ti-bands,  with  modifications  to  the  band  narar.eters 
that  can  be  measured  quantitatively. 

At  constant  “hu,  the  data  in  Fig.  2  and  in  "fan  charts"  constructed  for 
other  .samples,  are  well  explained  by  a  simple  two-band  model  .ipprux ima:  ion. 
Nevertheless,  t!ic  mass  parametei!.  determined  from  such  a  model  exhibi  t  :  i  i-'e 
dependence  on  liu,  slrcc  the  Land.iu  level  transitions  In  Fig.  2  show  de’w.r- 
tures  from  e  linear  dependence  between 'i*.!  and  resonant  magnetic  field.  I'lerc 
fore  the  comp.irison  between  masses  •'or  the  Interrnl  nte<f  and  pure  »■  r.iphl  te 
eemplrs  is  made  al  .t  common  reiorence  nulMr,  taken  at  the  K-noinf  '  lu’  oi'ci  , 
and  found  by  t.iklng  liw  *  U.  !!i«'  ^e^.>tlts  fur  the  vaJeiu*.'  and  lmiuJih  i  on  i.r  i,' 


Mognetic  field  (Tesla) 
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nnsses  n'  and  thus  ohCcninod  nro  rjiven  in  Table  1  for  a  number  of  dif- 
C  V 

ferent  samples.  We  note  that  in  ail  cases,  the  mass  parameters  in  the  inter- 
Tablc  1  Results  of  analysis  of  magnetoreflection  experiments 


Intercalant 

Stage 

* 

m 

c 

ft 

m 

v 

(Yq/Yj) 

(HOPG) 

oa 

0.056 

0.084 

25.1 

AlCl, 

8 

0.056 

0.076 

26.0 

6 

0.054 

0.076 

26.6 

FeCl, 

7 

0.055 

0.079 

26.1 

J 

5 

0.054 

0.075 

26.6 

Rb 

6 

0.045 

0.065 

31.7 

calation  compounds  differ  from  the  masses  in  graphite  by  only  a  few  percent, 
with  significantly  larger  changes  in  m“  and  m^J  occurring  in  the  donor  Rb  com¬ 
pound  than  in  acceptors  of  comparable  stage. 

Analysis  of  the  magnetoreflection  spectra  in  Fig.  2,  and  in  "fan  charts" 
obtained  for  other  samples  with  n^4  show  that  these  spectra  can  be  explained 
by  dispersion  relations  which  have  the  same  form  as  the  three-dimensional 
Slonczewski-Weiss-McClure  (SWIcC)  band  model  for  graphite,  but  with  modified 
values  of  the  band  parameter^.  Included  in  Table  1  are  values  for  the  SV.llcC 
band  parameter  combination  Yq/Yj^.  which  is  directly  related  to  the  reduced 
effective  mass  m’'  [4]. 

Magnetoreflection  resonances  such  as  in  Fig.  1  are  not  observed  below  a 
certain  photon  energy,  the  cutoff  energy  which  depends  on  both  stage 
and  intercalate  species.  The  cutoff  of  K-point  interband  transitions  bolow 

is  Interpreted  as  due  to  the  introduction  of  carriers  by  the  intercalant, 
resulting  in  a  shift  in  Fermi  level  Ep.  Because  of  the  Pauli  exclusion  prin¬ 
ciple,  interband  Landau  level  transitions  are  made  from  occupied  valence  to 
unoccupied  conduction  states.  The  lowering  of  Ep  caused  by  the  introduction 
of  holes  in  acceptor-type  compounds  leads  to  a  cutoff  of  Landau  level  tran¬ 
sitions  as  Ep  drops  below  the  extremum  of  the  magnetic  subbands  for  the  init¬ 
ial  state.  With  increasing  intercalate  concentration,  Ep  moves  to  lower  en¬ 
ergies,  therefore  increasing  the  magnitude  of  In  donor-type  compounds 

the  introduction  of  electrons  causes  Ep  to  rise,  leading  to  the  cutoff  of 
Interband  transitions  as  Ep  rises  above  the  extremum  of  the  magnetic  subbands 
for  the  final  state.  Analysis  of  the  measured  cutoff  photon  energies  yields 
the  shift  in  the  Fermi  level  relative  to  the  K-peint  band  edge  jEp-E^'^i.  . 
Results  for  several  intcrcalants  and  stages  are  presented  in  Fig.  3,  vMerc 
the  +  sign  in  Ep-Ej  is  for  donors  and  the  -  sign  for  acceptors.  Of  sig¬ 
nificance  is  the  much  larger  shift  in  the  Fermi  level  in  donor  compounds 
relative  to  acceptor  compounds  of  similar  stage,  indicating  n  signif ic.’ntlv 
larger  charge  transfer  in  the  case  of  donor  compounds.  To  carry  out  .<  oom- 
plcto  analysis  of  the  mngnetorcflcction  spectra,  it  is  necessary  to  u.-.e  Lan¬ 
dau  levels  for  intercalated  graphite,  and  this  is  discussed  below  In  Section 
5. 
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Fig.  3  Dependence  on  reciprocal  stage  (1/n)  of 
the  Fermi  level  measured  relative  to  the  K-noint 
band  edge  for  several  intercalants.  The  -r  sign 
for  (Ep-E^  applies  to  the  donor  compounds 
and  the  -  ^ign  to  the  acceptors.  For  all  in- 
tercalants,  the  graphite  value  is  obtained  in 
the  limit  (l/n)->0  and  is  denoted  by  open  squares. 
The  open  triangle  7  is  a  Shubnikov-de  Haas  re¬ 
sult  from  the  present  work  and  the  closed  tri¬ 
angle  T  is  from  magnetic  susceptibility  results 
of  DISALVO  ET  AL  [12]. 


3.  Shubnikov-de  Haas  Experiments 


Having  established  that  the  electronic  levels  near  the  Fermi  level  are  closely 
Identified  with  the  graphite  ir-bands,  it  is  of  interest  to  explore  the  effect 
of  intercalation  on  the  Fermi  surface.  In  this  connection  Shubnikov-c'e  Haas 
(SdH)  and  other  quantum  transport  experiments  have  been  reported  on  a  wide 
variety  of  graphite  intercalation  compounds  [5-8].  The  multitude  of  experi¬ 
mental  results  suggest  that  with  careful  s^lranle  preparation  and  character¬ 
ization,  a  complete  set  of  frequencies  associated  with  each  stage  o’*  a  aiven 
intercalation  system  can  be  obtained.  A  number  of  careful  experiments  have 
already  been  reported.  Perhaps  the  most  complete  set  of  data  has  been  ob¬ 
tained  for  the  graphite-K  system  where  SHEMATSU  ET  AL  [6]  give  results  for  stages 
n  ■  1,3,4  and  we  for  n  =  4,5,8.  A.lso,  results  for  graphite-Rb  compounds 

have  been  previously  reported,  showrrig  that  the  Fermi  surfaces  are  not  only 
stage-dependent,  but  also  vary  from  one  alkali  metal  intercalant  to  another. 


To  obtain  quantitative  information  on  the  stage  dependence  of  these  Fermi 
surfaces,  it  is  important  to  work  with  well-staged  and  characterized  sam.nl os. 
Since  the  Fermi  surfaces  for  the  alkali  metal  compounds  with  the  intercal- 
ants  K  and  Rb  are  stage-dependent,  staging  fidelity  is  essential  for  obtain¬ 
ing  reproducible  experimental  results  for  different  samples  of  the  same 
stage  index.  Because  of  the  instability  of  alkali-metal  samples  in  the 
presence  of  air  and  moisture,  the  samples  are  encapsulated  in  ampoules  and 
sample  handling  is  done  in  an  Argon-filled  dry  box  ('v^l  ppm  oxygen  coi'.uent,' . 

The  stage  of  the  samples  was  determined  using  (002.)  x-rny  diffraction  oru- 
files  both  before  and  after  the  SdU  experiments,  confirming  tlnit  t!ie  sa.:r'lcs 
were  single-staged  and  that  no  desorption  had  occurred  during  the  measure  lents . 

The  four-point  method  was  used  to  study  the  in-plane  transverse  mnvneto- 
resistence  In  the  temperature  range  1.4sTs4.2K  and  in  magnetic  fieiU.s  up  to 
15  Telsa.  The  leads,  were  attacheil  to  the  sample  using  cuiniucting  cnixy.  ;'ii 
sample  was  then  insirted  in  .i  Iielinm-Mlleu  ampoule  uu!  styc.isL  w.i  .  i>>  i  tu 
seal  the  ampoule.  The  angular  licnendeiire  of  tlie  S«jl!  ii  lat  j  uns  I  u  ■. 

be  me.j.siired  by  rotating  tl>e  namiile  around  the  ilircction  of  lac  curr.  c. 
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sm;!>  thnt:  I  for  all  angles,  anc!  lienee  transverse  ma.gnetoresistenee  was 
alwiiys  measured.  Data  acquisition  vias  by  computer,  and  the  data  were  mani¬ 
pulated  to  obtain  a  Fourier  power  spectrum  of  resistance  vs.  1/H,  thereby 
yielding  the  frequencies  of  SdH  oscillations.  These  frequencies  are  related 
to  the  extremal  cross  sections  of  the  Fermi  surface.  The  Fourier  power  soec- 
tra  for  the  Sdll  frequencies  for  stages  4,5,8  graphite-potassium  are  prcsonteci 
in  Fig.  4. 

During  the  course  of  these  studies,  a  great  deal  of  care  was  given  to  en¬ 
sure  tlie  fidelity  and  reproducibility  of  the  data.  For  this  reason,  the  ex¬ 
periment  Vvias  performed  on  five  potassium  stage  5  samples  and  on  one  of  these, 
the  measurements  were  repeated  several  weeks  later  after  attaching  new  leads. 
liThile  x-ray  profiles  after  each  experimental  step  confirmed  that  the  samples 
retained  their  single  stage  identity,  the  magnetoresistance  oscillations  and  - 
the  relevant  Fourier  power  spectra  revealed  nearly  identical  traces  for  these 
stage  5  samples.  Measurements  done  on  potassium  samples  with  different  stages 
or  on  rubidium  samples,  however,  showed  distinctly  different  SdH  frequencies. 

In  all  cases  the  SdH  frequencies  were  very  different  from  those  found  in 
pristine  graphite.  Hence  our  conclusion  is  that  there  is  a  unique  set  of 
SdH  frequencies  and  Fermi  surfaces  associated  with  each  stage  and  donor  in- 
tercalant  (K,Rb).  SUEMiTSU  ET  AL  [6]  have  also  shoiim  stage  dependent  de  Haas-van 
Alphen  frequencies  in  stages  3  and  4  graphite-potassium  compounds. 


Fig.  4  Shubnikov-de  Haas  Fourier 
Transform  Power  Spectra  for  stages 
4,  5  and  8  graphite-K.  These 
power  spectra  were  obtained  by  a 
Fourier  transform  of  an  experi¬ 
mental  resistivity  1/H  trace 
for  magnetic  fields  H  <  15  Tesla. 
The  peaks  in  the  power  spectra 
correspond  to  SdH  frequencies, 
which  are  given  in  Tesla  and  the 
same  scale  is  used  for  each  stage. 


Poorly-staged  samples  also  give  distinct  SdH  frequencies.  In  fact  more 
frequencies  are  found  in  mixed  than  in  single-stag  d  samples.  Sainnlcs  .-hew¬ 
ing  the  greatest  stage  fidelity  also  show  the  simplest  SdH  spectra.  I  iu  c  - 
fect  of  mixed  staging  introduces  additional  frequenc.'os  associated  with  f  ■ 
minority  phase  and  shifts  the  Fermi  level  to  some  intermediate  value  bctwifp 
chose  corresponding  to  t!\e  pure  stage  values.  This  shift  in  Ferni  leve'  }h-. 
most  sensitively  observed  in  the  very  low  frequency  cross-sections  where  a 
small  change  in  Fermi  level  can  result  in  large  shifts  in  frequency.  There 
is  also  some  indication  that  different  in-plane  intercalate  densities  will 
yield  shifts  in  Fermi  level.  In  fact,  a  detailed  Interpretation  of  the  St’!! 
spectra  is  expected  to  provide  a  critical  test  for  any  model  of  tlie  elee- 
Cronic  structure  for  tliese  matt  i  lals:. 
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Measurements  ol  tlie  iiiit;ul.:r  dependence  of  the  Shiihnikov  de  Hans  f*-en"er>.- 
cies  (shown  in  Fig. 5)  suggest  that  the  Fermi  surface  is  segmented  into 
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Fig.  5  Angular  dependence  Q(C)  /?.(C)  c~ 
several  dominant  SdH  frequencies  wlierc- 
corresponds  to  H;,  c-axis  and  th.e  in¬ 
dicated  frequency  values  are  for  j=0. 

The  290,  191  and  24  T  frequencies  ere 
not  observed  above  '^30®,  n.45®  and  ■'■>15°  , 
respectively.  The  angular  dependence  of 
the  290  T  and  191  T  frequencies  is  witl;- 
in  experimental  error  given  by  a  cos  6 
dependence  (solid  line),  characteristic 
of  a  cylindrical  Fermi  surface,  x.’hile  Che 
24  T  frequency  exhibits  departures  from 
the  cos  6  dependence  at  large  angles. 


cylindrical  rings.  Deviations  from  two-dimensional  behavior  is  measured  by 
small  departures  from  the  cylindrical  angular  dependence  and  by  the  obser¬ 
vation  of  doublet  SdH  frequencies  (Fig. 4)  which  indicate  different  cross- 
sectional  areas  at  the  K  and  H  points  in  the  Brillouin  zone. 

4.  Application  of  Enerev  Band  Model 


To  Interpret  these  Shubnikov-de  Haas  results  a  band  model  appropriate  to  a 
range  of  intercalate  stages  is  necessary.  Such  a  model  for  the  electronic 
dispersion  relations  F('k)  for  graphite  intercalation  compounds  is  clready 
available  in  terms  of  a  Hamiltonian  based  on  the  S'.^’McC  three-dimensio:-..i’ 
model  for  the  graphite  rr-bands  [9].  In  this  model  Che  c-axis  sunerlattite 
periodicity  is  included  through  a  k^-axis  zone-folding  of  the  energy  levels. 

A  unitary  transformation  then  transforms  the  matrix  Hamiltonian  into  n  layer 
representation.  For  a  stage  n  compound,  n  graphite  layers  are  retained  and 
the  (n+l)st  layer  is  replaced  as  a  first  approximation  by  an  empty  interca¬ 
late  layer  and  for  more  quantitative  results  by  an  intercalate  layer  v;bich 
interacts  with  the  adjacent  graphite  bounding  layers.  In  Che  "emuty  inter¬ 
calate  layer"  model,  the  matrix  Hamiltonian  depends  only  on  the  graphite  band 
parameters,  and  these  are  already  known  from  previous  experiments  on  pristine 
graphite.  The  raagnetoreflection  experiment  provides  revised  values  for  some 
of  these  band  parameters,  as  modified  by  the  intercalation  process.  Intcr- 
calant-specif ic  interactions  between  the  intercalant  and  the  graphite  bound¬ 
ing  layer  can  then  be  introduced  to  obtain  the  final  dispersion  relations. 

This  model  for  the  electronic  dispersion  relations  has  been  applied  to  the 
interpretation  of  th:  observed  SdH  f renquencies .  Since  the  most  conp.lece  sc*. 
of  experimental  data  is  available  for  stage  5  graphite-K,  the  aoplicat ion  o'’ 
the  energy  band  model  is  made  for  this  case.  Thus  in  Fig.  6  are  sln'v.n  '.  '.e 
electronic  dispersion  relations  for  the  five  valence  and  five  conduction  - 
bands  appropriate  to  a  stage  5  intercalation  compound.  The  emp Ir.ical  position, 
of  the  Fermi  level  Ep  shov.'n  in  Fig.  6  is  adjusted  to  yield  the  best  fit  ’*0  the 
experimental  SdH  frequencies,  in  this  case  corresponding  to  four  partial'.v 
occupied  conduction  bands.  The  Fermi  surface  parameters  obtained  from  t’lis 
model  are  listed  in  Tab,le  2,  including  calculated  anu  observed  SdH  friq'anc- 
les  :’(0)  for  111'  c-nxis,  the  cyclotron  mass  at  lip,  ilic  trigonal  waruin/, 

anisotropy  listed  as  Uj/k2,  -aid  tlie  i-loctron  density  iij  lor  eaca  c.arrii.r 
pocket.  The  values  for  k^^/k^  in  Table  ?.  suggest  that  trigonal  warping  ;s 
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Fig.  6  Electronic  band  structure  in  the 
vicinity  of  the  HK  axis  near  the  Fermi  level 
for  a  stage  5  "empty  Intercalate  layer"  model 
for  a  graphite  intercalation  compound.  The 
Fermi  level  is  determined  to  fit  the  ob¬ 
served  SdU  frequencies.  For  the  indicated 
Fermi  level  four  of  the  conduction  bands  are 
partially  occupied  and  give  rise  to  Fermi 
surfaces  and  carrier  pockets. 


Table  2  Fermi  surface  parameters  associated  with  stage  5  graphite-K 


Fermi  Surface 
Parameters 

*^1 

K-point  Band  Designations 

•  m 

0.143 

0.115 

0.0828 

0.0511 

n^(xlO^°cra“^) 

2.36 

1.76 

0.975 

0.138 

Calculated  SdH 

Frequencies  n(0) 

401 

300 

163 

26.7 

Observed  {2(0) 

453 

290 

191 

.24  1 

(Tesla) 

430 

267 

152 

18 

243 

135 

Trigonal  Harping 

Anisotrophy  kj^/k2 

0.68 

0.77 

0.91 

1.00  1 

important  for  the  larger  cross-sectional  areas  with  heavier  masses,  whereas 
the  smaller  light  mass  cross-sections  are  circular.  The  generally  good  agree¬ 
ment  of  the  observed  IdH  frequencies  with  the  empty  intercalate  Icye^-  modal 
of  Fig. 6  suggests  that  the  effect  of  intercalant-graphite  bounding  layer  in- 
tcractluns  can  be  treated  as  a  perturbation  and  evaluated  by  fitting  the  model 
quantitatively  to  the  observed  Sdil  frequencies.  It  should  be  poted  chat  th.- 
volumes  of  the  carrier  pockets  Cor  the  four  occupied  bands  in  Fig.  6  corres¬ 
pond  to  a  charge  transfer  of  ’v>0.3  electrons  per  intcrcalant  into  these  carrier 
pockets,  suggesting  that  other  carrier  pockets  could  be  present  elsewhere  in 
Che  Briilouin  zone. 

5.  T..indau  Level  Calculation 


The  empty  iiitercnlatc  layer  model  used  Co  calculate  the  Fermi  surface  can 
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also  be  applied  to  compute  the  Landau  levels.  The  results  of  such  a  calcu¬ 
lation  are  shovm  in  Fig.  7  for  a  stage  3  intercalation  compound.  For  a 
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Fig.  7  Magnetic  field  dependence  of  the 
four  lowest  quantum  number  Landau  levels 
at  kj,"®  ^  stage  3  graphite  intercala¬ 

tion  compound,  calculated  on  the  basis  of 
an '(empty  intercalate  layer"model.  The 
middle  set  of  levels  correspond  to  the  E. 
•levels  of  graphite  while  the  upper  and  lower 
sets  correspond  respectively  to  the  and 
E2  levels. 


stage  3  compound,  model  calculations  yield  three  conduction  bands  and  three 
valence  bands,  such  that  two  of  the  three  conduction  and  valence  bands  lie 
close  to  each  other  in  the  vicinity  of  the  HK  axis,  leaving  one  higher-lying 
conduction  band  and  one  lower- lying  valence  band.  (See  Fig. 6  where  the  levels 
for  a  stage  5  compound  are  shorn.)  The  near-degeneracy  of  the  four  bands  Jor 
the  case  of  stage  3  results  in  magnetic  level  crossings  as  shown  in  Fig. 7. 

In  this  calculation  the  trigonal  warping  is  also  neglected  only 

linear  k  terms  in  the  off-diagonal  magnetic  Hamiltonian  are  retained.  Under 
these  assumptions  the  Landau  levels  can  be  calculated  exactly  and  the  usual 
optical  selection  rules  apply  (An=±l).  The  central  set  of  levels  in  Fig.  7 
correspond  to  the  K  and  H  point  transitions  between  the  pristine  graphite 
Ej  bands,  including  levels  arising  from  k  -axis  zone  folding.  On  the  ocher 
hand,  the  upper  and  lower  Landau  levels  of  Fig.  7  are  associated  with  the 
E]^  and  E2  graphite  bands.  Landau  level  transitions  involving  these  bands 
are  analogous  to  the  H-polnt  transitions  in  pristine  graphite. 


Assuming  that  the  Fermi  level  in  a  stage  3  donor  compound  lies  between  the 
central  and  upper  group  of  magnetic  energy  levels  (see  Fig. 7),  additional 
sets  of  Landau  level  transitions  are  predicted,  corresponding  to  the  excita¬ 
tion  of  an  electron  from  the  central  set  of  levels  to  the  upper  set.  Such 
excitations  require  a  photon  energy  of  4.0.6  eV  which  is  in  excess  of  the  pho¬ 
ton  energies  employed  to  date  in  the  magnetoreflection  experiments  on  inter¬ 
calated  graphite.  The  Landau  level  calculation  itself  is  vital  for  the  de¬ 
tailed  analysis  of  the  magnetoreflection  experiment  with  regard  to  identi¬ 
fication  of  the  observed  resonant  structures,  selection  rules  for  Landau 
level  transitions  and  lineshape  analysis  for  determination  of  the  resonant 
point  within  the  experimental  linewidth.  To  carry  out  such  a  detailed  an¬ 
alysis  it  is  necessary  to  include  explicitly  the  Interaction  between  the 
intercalant  and  the  graphite  bounding  layers.  In  turn,  such  an  analysis 
provides  a  sensitive  method  for  the  experimental  determination  of  the  mag¬ 
nitude  of  these  interaction  parameters. 

The  Shubnlkov-de  Haas  oscillations  in  the  conductivity  are  also  related 
to  the  Landau  levels  shown  in  Pig.  7.  However,  for  those  cross  sections 
with  fast  frequencies  where  trigonal  warping  (sec  Table  2)  is  important,  it 
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is  necessary  to  also  include  explicitly  the  trigonal  \jarning  parnc^cCor 
and  higher  terms  in  k  In  order  to  make  a  quantitative  fit  to  be  experimental 
SdH  frequencies.  A  systematic  fit  of  both  Shubnlkov-de  Haas  and  magneto- 
reflection  data  for  a  single  series  of  compounds  is  now  underway  and  is  ex¬ 
pected  to  result  in  precise  determination  of  the  interaction  between  the  in- 
tercalant  and  the  graphite  bounding  layers. 
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